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ABSTRACT 

We estimate the relative contributions of the supermassive black hole (SMBH) accretion disk, corona, 
and obscuring torus to the bolometric luminosity of Seyfert galaxies, using Spitzer mid- infrared (MIR) 
observations of a complete sample of 68 nearby active galactic nuclei (AGNs) from the INTEGRAL 
all-sky hard X-ray (HX) survey. This is the first HX-selected (above 15 keV) sample of AGNs with 
complementary high angular resolution, high signal to noise, MIR data. Correcting for the host galaxy 



contribution, we find a correlation between HX and MIR luminosities: L 



i5 M mOcL^ 4±u ub . Assuming 
that the observed MIR emission is radiation from an accretion disk reprocessed in a surrounding dusty 
torus that subtends a solid angle decreasing with increasing luminosity (as inferred from the declining 
fraction of obscured AGNs), the intrinsic disk luminosity, ^Disk, is approximately proportional to the 
luminosity of the corona in the 2-300 keV energy band, ^corona, with the lvDisk/-kcorona ratio varying 
by a factor of 2.1 around a mean value of 1.6. This ratio is a factor of ^2 smaller than for typical 
quasars producing the cosmic X-ray background (CXB). Therefore, over three orders of magnitude in 
luminosity, HX radiation carries a large, and roughly comparable, fraction of the bolometric output 
of AGNs. We estimate the cumulative bolometric luminosity density of local AGNs at ~(1 — 3) x 
10 40 erg s _1 Mpc -3 . Finally, the Compton temperature ranges between kT c « 2 and « 6 keV for 
nearby AGNs, compared to kT c w 2 keV for typical quasars, confirming that radiative heating of 
interstellar gas can play an important role in regulating SMBH growth. 

Keywords: galaxies: active - galaxies: Seyfert - infrared: galaxies - X-rays: galaxies 



. 1. INTRODUCTION 

Active galactic nuclei (AGNs) are extremely powerful 
(N. sources of electromagnetic radiation over many decades 
in frequency from radiowaves to gamma-rays. Accord- 
. ing to the commonly accepted scenario, an AGN shines 
due to accretion of gas onto a supermassive black hole 
(SMBH) residing in a galactic nucleus. 
In Seyfert galaxies and quasars, most of the luminos- 
■ ity is emitted in the form of ultraviolet (UV) radiation 
' generated in a geometrically thin, optically thick accre- 
tion disk (|Shakura fc Sunvaevlll973l ). giving rise to a "big 
• i— i . blue bump " (BBB) in the spectral energy distribution 
^ ■ (SED, e.g. JMalkan fc Sargentlll982l h Additional, higher 
energy radiation is generated in a hot coro na of the ac- 
?3 ■ cretion disk fe.g.. lHaardt fc Maraschilll993f) and possibly 
also in collimated outflows (jets) of relativistic plasma, 
producing a hard X-ray (HX) peak in the SED. The in- 
tegrated (and redshifted) HX emission of all AGNs in 
the observable Universe makes up the bulk of the cosmic 
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X-ray background. There is als o a third, mid- infrared 
(MIR) peak in AGN SEDs (e.g.. lBarvaimsl[T987t ). which 
arises from reprocessing of a significant fraction of the 
disk's and some of the coronal radiation in a torus of 
molecular gas and dust surrounding the inner accretion 
flow. In fact, only in unobscured or "type 1" AGNs can 
all three spectral components, the HX bump, the BBB, 
and the MI R bump, be obser ved. According to the uni- 
fied model (jAntonucci 111993ft . these are objects viewed 
through the funnel of the dusty torus. In contrast, only 
the HX and MIR components are visible in the SEDs of 
obscured or "type 2" AGNs, because the torus is opaque 
to UV emission from the accretion disk but transparent 
to coronal radiation at energies above ~15 keV (except in 
Compton-thick sources) and to its own infrared emission 
(at least at wavelengths > 20 /im). All other emission 
components, including broad- and narrow-line emission 
and non-thermal radio and gamma-ray radiation are usu- 
ally not significant as regards their contribution to the 
angular-integrated bolometric luminosity of AGNs; these 
components will therefore not be discussed below. 

To understand how electromagentic radiation is emit- 
ted and reprocessed during accretion of matter onto 
SMBHs, it is crucial to explore i) in what proportion the 
AGN luminosity is shared between the accretion disk and 
its corona, ii) what fraction of the bolometric luminosity 
is reprocessed in the torus, and iii) how these proper- 
ties depend on black hole mass and accretion rate. One 
also needs such information to study the role of AGN 
feedback in regulating SMBH growth and galactic evolu- 
tion. One of the proposed feedback mechanisms is pho- 
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toionization and Co mpton heating of interstellar gas by 
AGN radiation (e.g.. lCiotti & Ostrikerl2"00lHProga et all 
I2008D . whose efficien cy cr i tically depends on the AGN 
SEP (jSazonov et al.1 [200l 120051 ) . Finally, information 
on AGN SEDs can be used to derive bolometric cor- 
rections required to reconstruct the cosmic history of 
SMBH accretion growth based on AGN statistics pro- 
vided by extragalactic surveys (e.g.. lMarconi et aT]|2004t 
iMerloni fc HeirI3l2008h . 

Among all types of AGNs, the SEDs of unobscured 
high-luminosity quasars have been studied most exten- 
sively (see, e.g.. [Elvis et al.lll994t iRichards et"aH 120061 : 
iShang et al.ll2011[ ). Their obscured counterparts - type 2 
quasars - have been explored to a much lesser degree, 
although recent surveys have beg un to find such ob- 
jects in significant numbers (e.g. , iPolletta et al.1 120061 : 
iHickox et all 120071 : lLanzuisi etall I2009Q . There is also 
much uncertainty with respect to the SEDs of Seyfert 
galaxies, which are typically less luminous than more 
distant quasars. The difficulty is that even in Seyfert Is, 
the accretion disk emission is usually contaminated by 
host galaxy stellar emission in visible bands and the BBB 
peaks in the observationall y difficult far-UV band (see , 
howe ver, iScott et all l2004t iVasudevan fc Fabianl 120071 
120091) . 

The goal of the present study is to systematically as- 
sess the relative contributions of the accretion disk, hot 
corona, and obscuring torus to the bolometric luminos- 
ity of local Seyfert galaxies. To this end, we i) cross- 
correlate the HX luminosities of nearby AGNs detected 
during the all-sky survey of the International Ga mma- 
Ray Laboratory (INTEGRAL, \Wmk\er et al. II2003D with 
the MIR luminosities of these objects meas ured by the 
Spitzer Space Telescope (jWerner et al.1120041 ). and ii) use 
the proportion of obscured to unobscured AGNs to esti- 
mate the opening angle of dusty tori as a function of lu- 
minosity. We then put our findings for nearby AGNs into 
the broader context of cosmic SMBH growth by making 
a comparison with distant quasars. 

Most previous relevant studies were based on AGN 
samples compiled in a fairly arbitrary manner from op- 
tical and/or soft X-ray (below 10 keV) catalogs (e.g ., 
iLutz et al] [2001 iHorst et all [200l iHonig et al] 1201(1 . 
In these energy bands, AGNs can easily be missed due 
to absorption, as powerful sources can become invisible 
when obscured by large amounts of dust and cold gas 
in the torus and/or host galaxy. Furthermore, as al- 
ready noted above, optical emission from relatively low- 
luminosity AGNs can be diluted against th e background 
of a luminous galaxy (see [Mushotzky 2004 for a detailed 
discussion of AGN selection at different wavelengths). 

The hard X-ray band, above ~ 15 keV, provides a 
census of AGNs that is far less biased with respect to 
the viewing orientation of the torus and is unbiased 
with respect to host galaxy p roperties. There have 
been a few previous attempts (jVasudevan et al] 120101 : 
iMullanev et al.|[20Tl|) of systematically studying the MIR 
properties of HX s elected AGNs using the Swift all-sky 
hard X-ray survey (|Tueller et al.l 120101 ) . However, these 
studies either used data from the IRAS all-sky photo- 
metric infrared survey, so that it was impossible to re- 
liably subtract the host galaxy contribution from the 
AGN emission, or used high angular resolution Spitzer 
data but only for statistically incomplete subsamples of 



Swift AGNs. Our INTEGRAL sample is the first statis- 
tically complete, HX selected sample of AGNs with com- 
plementary high angular resolution, high signal to noise, 
MIR data. The extensive Spitzer coverage (3.6-38 /im) 
available for the entire INTEGRAL sample makes this a 
unique data set for studying SEDs of AGNs in the local 
Universe. 

2. INTEGRAL AGN SAMPLE 

Our study is based o n the complete sam ple of AGNs 
(|Krivonos et al. 1 120071 : iSazonov etall I2007D detected in 
the 17-60 keV energy b and by the IBIS/ISGRI detec- 
tor (jUbertini et al.ll2003h aboard INTEGRAL during the 
first three and a half years of the mission, from 2002 Oc- 
tober until 2006 June. These observations compose a 
serendipitous all-sky HX survey with the flux limit vary- 
ing by a factor of a few over the sky. We have excluded 
from the present analysis blazars (flat spectrum radio 
quasars and BL Lac objects), a relatively rare subclass 
of AGNs whose observed emission is believed to be dom- 
inated by a narrow, strongly collimated component. We 
have also excluded AGNs located in the "zone of avoid- 
ance" near the Galactic plane (\b\ < 5°) because there 
remain unidentified INTEGRAL sources in this region 
of the sky while we wish our sample to be nearly 100% 
complete to minimize selection effects. 

The resulting set comprises 68 AGNs (Table [lj. In 
the first seven columns of Table Q] we have collected in- 
formation on optical/radio AGN types, distances, HX 
fluxes and luminosities, and X-ray absorption column 
densities (-/Vh)- These data are mostly adopted from 
the original INTEG RAL catalog (jKrivonos et al. 1 120071 : 
ISazonov et al.ll2007l ) although some updates take into ac- 
count follow-up observations carried out since publica- 
tion of the catalog. In particular, thanks to recent X- 
ray observations by Chandra, Swift, and XMM-Newton, 
all of the previously missing TVh values have now been 
estimated. All the reported absorption columns may be 
considered reliable because they are based on high signal- 
to-noise X-ray spectroscopic data. We do not quote the 
uncertainties associated with the iVn values, because the 
information on absorption columns has been compiled 
from various sources and in most cases the actual un- 
certainty is likely dominated by systematic effects asso- 
ciated with the particular spectral modeling procedure 
used. In fact, multiple measurements taken for some 
AGNs at different times and/or by different instruments 
sometimes yield values that differ from each other 
by more than their reported uncertainties. We estimate 
that the total uncertainties associated with A^h columns 
for our sources are typically smaller than 30% and do not 
affect the present study in any significant way. 

All but one of our AGNs are located at low redshift 
(z < 0.1, the most distant one, IGR J09446-2636, be- 
ing at z — 0.14). For 18 nearby (closer than ~ 40 Mpc) 
Seyfer t galaxies we have ado pted di s tance estimates from 
either EJEitaD (pOOl ) or lTulTvl (fl988h : otherwise lu- 
minosity distances have been calculated from the spec- 
troscopic redshifts assuming a cosmology with Q m = 0.3, 
Q A = 0.7, and H = 72 km s" 1 Mpc" 1 . 

The INTEGRAL AGN sample has HX (17-60 keV) lu- 
minosities ranging over almost five orders of magnitude, 
from 4 x 10 40 (NGC 4395) to 2 x 10 45 (IGR J09446-2636) 
erg s . Therefore, this is a representative, hard X-ray 
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Table 1 

INTEGRAL-Spitzer AGN sample 



Name 



AGN a 
Class 



Ref 



D, 
Mpc 



Ref 



10 22 cr 



Ref 



17-60 keV 

TT 



15 /an 



Flux, 10 
erg s — 1 cm" 



logi, 



erg s 



Jy 



\o%vL v , Fagn 
erg s — 1 



Clean sample (AGN dominated infrared sources) 
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selected sample of nearby AGNs, mostly Seyfert galax- 
ies, although our ~ 10 most luminous objects may be 
better referred to as nearby quasars, because their HX 
luminosities exceed 10 44 erg s -1 . 

One special object in the sample is the Seyfert 
1.8 galaxy NGC 4395, a famous low-luminosity AGN 
(LLAGN) sometimes referred to as a "dwarf Seyfert nu- 
cleus" . It appears to be a quite atypical Seyfert galaxy in 
terms of its black hol e mass, luminosity , and variability 
properties (see, e.g.. iMoran et al.l 119991 : iPeterson et al.l 
120051 : IVaughan et al.ll2005D . It is therefore possible that 
the properties of its dusty torus (if there is any) are also 
different from typical Seyfert galaxies and quasars. We 
thus treat this object separately from the rest of the sam- 
ple in performing the HX-MIR cross-correlation analysis 

As noted above, the HX selection (17-60 keV) implies 
that there is almost no bias from absorption. Our AGN 
sample is not sensitive to photoabsorption (and we have 
thus not corrected the measured HX fluxes for line-of- 
sight absorption) as long as the column density of the gas 
is less than a few 10 24 cm~ 2 or equivalently the Thom- 
son optical depth is less than a few; at even larger col- 
umn densities, the flux from a source drops considerably 
at all X-ray energies. The Seyfert 2 galaxy NGC 1068 
distinguishes itself from the rest of the sample because 
it is the only significantly Compto n-th ick AGN (hav- 
ing N H > 10 25 cm" 2 , iMatt et al.l 120001 ) . We therefore 
exclude NGC 1068 from our baseline HX-MIR cross- 
correlation analysis but discuss its properties in compar- 
ison with Compton-thin sources ([JSJ 

3. SPITZER OBSERVATIONS AND DATA REDUCTION 

More than half of the INTEGRAL sample consists of 
well-known Seyfert galaxies, many of which have been 
targets of observational campaigns with Spitzer. For the 
remaining part, largely represented by AGNs discovered 
by INTEGRAL, we carried out short Spitzer observations 
(Program ID 50763) consisting of 3.6-8 /im imaging with 
IRAC and low-resolution MIR spectroscopy with IRS; 
in addition, far-infrared photometry was performed with 
MIPS for a subset of objects. Our proprietary and pub- 
licly available archival data together provide complete 
coverage of the INTEGRAL sample with Spitzer at 3.6- 
38 /im. 

All of the sources in our hard X-ray selected sample 
would have been robustly detected by Spitzer even if they 
had been 1-3 orders of magnitude fainter. Hence, our 
sample is not limited by MIR flux and any correlations 
derived between the HX and MIR luminosities can be 
considered representative of the local AGN population 
without significant bias. 

3.1. IRS 

We used the InfraRed Spectrograph fIRS. lHouck et al.l 
2004) on Spitzer to obtain low-resolution spectra of our 
objects. Our program's observations, for a total of 30 
AGNs, were done in spectral mapping mode using the 
Short-Low (SL) and Long-Low (LL) IRS modules. Each 
of these modules has first- and second-order sub-slits 
(SL1, SL2, LL1, and LL2) with widths of 3.7, 3.6, 10.7, 
and 10.5 arcsec, respectively. We used SL1, LL1, and 
LL2 for the entire sample and SL2 for a subset of objects. 



The resulting spectra thus cover a range from either 5.2 
or 7.5 fim up to 38 /im. Observations with SL1 and SL2 
consisted of one cycle of 6 pointings with a ramp duration 
of 6 s with two 19" steps in the slit direction and three 
1.8" steps in the dispersion direction. For LL1 and LL2, 
one cycle of 3 pointings with a ramp duration of 6 s and 
a step size of 42" in the slit direction was implemented. 

For those sources that were not covered by our Spitzer 
program, we used archival low-resolution IRS data: map- 
ping mode observations for 14 AGNs and staring mode 
observations for another 23 AGNs. Almost all of the 
archival mapping mode observations have the following 
setup: one cycle of 13 pointings with a step size of 1.8" 
in the dispersion direction for SL1 and SL2 and 5 point- 
ings with a step size of 4.85" in the dispersion direction 
for LL1 and LL2, the ramp duration being 6 s. For the 
archival staring mode observations, the ramp times and 
numbers of cycles vary from one object to another. 

In the analysis of mapping mode observations, we used 
basic calibrated data (BCD), extracted from the Spitzer 
Science Center pipeline (versions S18.0.1 and newer for 
our program's observations and versions S15.3.0 and 
newer for the archival observations). For each order of a 
given IRS module and for each position of a given source 
within the slit, we first produced a background image. 
For our program's observations, this was done by aver- 
aging over 2D spectra obtained in significantly (~ 19" for 
SL and ~ 42" for LL) off-source positions, whereas for 
the archival observations, a similar averaging was done 
over 2D spectra obtained in the other order of a given 
IRS module. The background image was then subtracted 
from the on-source 2D spectrum. Then, a ID spectum 
of the source was obtained using the Spitzer IRS Cus- 
tom Extraction software (SPICE) by applying the reg- 
ular extraction algorithm, which uses an aperture that 
gradually increases with wavelength in accordance with 
the telescope's point spread function (e.g., the aperture 
width is 7.2" at 6 /j,m and 36.6" at 27 /im). Finally, an 
averaging over a set of ID spectra extracted in different 
source positions within the slit was done. The staring 
mode spectra were obtained simply by averaging over 
nod-subtracted post-BCD spectra. These were derived 
with the same regular extraction algorithm as was used 
in our analysis of mapping mode observations. 

Although the data reduction procedure described 
above is fairly simplistic and does not fully exploit the 
potential of mapping mode observations (which are avail- 
able for nearly two thirds of our sample) , it is adequate 
for the purposes of our HX-MIR cross-correlation study, 
as confirmed by a comparison with an alternative, more 
detailed analysis of IRS data for a subset of INTEGRAL 
sources (see §4.4|) . 

The Compton-thick Seyfert 2 galaxy NGC 1068 is an 
extremely bright (~ 15 Jy at 15 /im) infrared source, 
which caused saturation of IRS. We therefore quote 
in Table Q] its MIR flux estimate based on a compi- 
lation of high angu lar resolution, infrared observations 
(jMason et al.ll2006t ). 

3.2. IRAC 

We used i mages obtained by the InfraRed Array Cam- 
era (IRAC, iFazio et al.1 [200l to determine source flux 
densities at 3.6, 4.5, 5.8, and 8.0 /im. This wavelength 
range partially overlaps with that covered by IRS spec- 
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troscopy, enabling direct comparison between the spec- 
troscopic and photometric results and providing an ex- 
tension of spectra to shorter wavelengths. 

Our IRAC program included 35 INTEGRAL sources. 
Observations were done in high dynamic range (HDR) 
mode. Specifically, a combination of a 0.6 s frame and a 
12 s frame was repeated in 9 random dithers for each 
source. For the rest of the sample, we made use of 
archival observations, for which the number and dura- 
tion of frames varied from source to source. Since most 
of our objects are bright infrared sources, our analysis 
was in most cases based on stacking the 0.6 s frames. 
For the 5 weakest sources (< 10 mJy at 3.6-8.0 /im), 
including LEDA 138501, ICR J09446-2636, ESO 263- 
G013, ICR J10386-4947, and NGC 4395, to improve the 
accuracy, we determined the fluxes by stacking the long 
(2, 12, or 30 s) frames. This was also done for 3C 111 
and Mrk 3 despite their relative brightness because the 
observations were not done in HDR mode. We verified 
that none of the sources was saturated. 

We analyzed post-BCD data using the standard point 
source extraction package APEX, part of the MOPEX 
software. Post-BCD images are adequate for this work 
because their main deficiency, poor artifact correction, is 
unimportant for 0.6-s frames. We estimated source fluxes 
by integrating the surface brightness in different aper- 
tures with radii between 2.4" to 12" and correcting for 
flux leakage outside the aperture under the assumption 
of a point-like source. Although this procedure is inaccu- 
rate for measuring fluxes of spatially extended sources, it 
is good enough to indicate the presence of extended host 
galaxy emission as a significant difference between fluxes 
measured in large and small apertures (see Appendix [A"]l . 

The existing IRAC photometric measurements for the 
Compton-thick Seyfert 2 galaxy NGC 1068 proved to be 
saturated and hence were not used. 

4. INFRARED SPECTRA: AGN MIR EMISSION 

Infrared emission from dust associated with star forma- 
tion in the host galaxy can provide a significant contri- 
bution to AGN MIR spectra. Therefore, to study torus 
emission, we need to estimate and subtract the star for- 
mation contribution from the Spitzer data. We observe 
clear signatures of star formation in many of our IRS 
spectra. These include polycyclic aromatic hydrocarbon 
(PAH) emission features and the vF v continuum rising 
toward the far-infrared. We therefore modelled the mea- 
sured spectra by a sum of starburst and AGN compo- 
nents, similarly to a number of previo us studies (e.g., 
iNetzer et alll2007t iMullanev et al.ll2011l ). 

Figure [T] shows examples of IRS spectra with negligi- 
ble, significant, and strong star formation contribution 
as deduced using the fitting procedure described below. 
At short wavelengths, we also show IRAC photometric 
fluxes measured in 2.4" and 12" apertures. 

4.1. Spectral decomposition using a starburst template 

We adopted the starburst template from lBrandl et all 
(20060, which is an average over low-resolution IRS 
spectra of a dozen nearby (D < 100 Mpc) starburst 
galaxies. This template is well suited for our analysis be- 

8 http://www.strw.leidenuniv.nl/~brandl/SB_template.html 



cause it was obtained by low-resolution IRS spectroscopy, 
similarly to the spectra studied here. 

We normalized starburst components in our objects 
based on the observed strength of PAH lines, which 
are believed to be a generic signature of star formation 
(|Roche et al. I [1991 . The presence of an AGN in a star 
forming galaxy may lead to a weakening of the PAH spec- 
tral features because PAH molecu les can be destroyed by 
hard AGN radiation (|Voitlll992[ ). However the impor- 
tance ofj^hisjsffejrt^^ .g.. | Smith et al.l 
12007b- lO'Dowd et all [200l iSales et al.ll201Q[ ). and we 
have assumed that the shape of the starburst spectral 
component is not affected by the presence of a central 
AGN. 

Our analysis consisted of the following steps. First, we 
fitted the spectra around (typically within ±0.6 /im of) 
the 6.2 /im and 11.3 /im PAH lines by a sum of a lin- 
ear continuum and a Gaussian. We then compared the 
derived PAH line fluxes with the corresponding values 
for the starburst template, which yielded two indepen- 
dent estimates of the amplitude of the star-formation 
component. The average of these two values was then 
adopted as the normalization of the starburst template. 
On average, the coefficients implied by the 6.2 /im and 
11.3 /im features proved to be consistent with each other, 
although there is ~ 40% scatter around the 1:1 ratio 
of the two coefficients. This indicates that there are 
~ 20% systematic uncertainties in the derived ampli- 
tudes of starburst components for our objects (this is- 
sue is further discussed in £ 14.41 below) . If there were 
no IRS data for the 6.2 /im feature (i.e., only first-order 
SL data at > 7.5 fim were available), we used the flux 
of the 11.3 /im line to normalize the starburst compo- 
nent. If the observed PAH features proved to be strong 
enough (depending on the source brightness, we required 
the PAH equivalent widths, EW, to be larger than 0.01- 
0.02 /im, as compared to EW=0.45 and 0.55 /im for the 
6.2 /im and 11.3 /im bands, respectively, in the starburst 
template) , we subtracted the estimated starburst contri- 
bution from the total spectrum to derive the AGN com- 
ponent. Otherwise, we considered star formation con- 
tamination insignificant and did not perform any sub- 
traction. 

We applied an additional correction to 13 spectra that 
contained detectable PAH features and exhibited a sig- 
nificant (>10-20%, depending on the source brightness) 
discontinuity near 14 /im, where the short-wavelength 
segment measured in the 3. 7"- wide SL1 slit connects to 
the long-wavelength segment measured in the 10.5"-wide 
LL2 slit. Such a "jump" is mostly likely caused by an 
extended source, i.e., it cannot be due to the AGN. In 
these objects, the long-wavelength (> 14 /im) part of the 
starburst component was rescaled to make the AGN com- 
ponent smooth across the SL-LL boundary. We did not 
make such a correction for NGC 4395, by far the weak- 
est infrared source in our sample (with an estimated flux 
density of 13 mJy at 15 /im), despite the apparent pres- 
ence of a significant SL-LL discontinuity in its spectrum, 
because of its low statistical quality (see Fig. [IJ. Fur- 
thermore, our comparison with available high-resolution 
spectroscopy for this object (see £14. 21 below) indicates 
that the SL slit was not positioned sufficiently accurately 
on the nucleus of NGC 4395, which might have caused 
an artificial discontinuity at 14 /im in the low-resolution 
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Figure 1. Examples of Spitzer spectra of INTEGRAL AGNs. In each panel, the solid line shows the low-resolution IRS spectrum, and 
the filled and open circles represent fluxes at 3.6, 4.5, 5.8, and 8.0 (im derived from IRAC images in 2.4" and 12" apertures, respectively 
(corrected for flux leakage outside the aperture assuming a point source). Also, the X-ray absorption columns arc indicated. Left column: 
Spectra that are clearly dominated by dust emission associated with the active nucleus, with a weak or absent starburst contribution. 
Middle column: Spectra showing a noticeable contribution of MIR emission from dust associated with star formation. The dotted line 
shows the starburst component, estimated by fitting the starburst template to the 6.2 and 11.3 (im PAH lines (indicated in the upper 
panel). The dashed line shows the AGN contribution, found as the difference between the total spectrum and the starburst component; 
it may still be contaminated at short wavelengths by stellar and accretion disk emission. Right column: Starburst dominated spectra, for 
which extraction of an AGN component is practically impossible. Note the much broader flux range covered by the spectrum of NGC 4945 
compared to the other sources. 
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IRS spectrum. 

In principle, we could also use another known strong 
PAH feature, at 7.7 /im, for estimating the contribu- 
tion of star formation. However, this band overlaps with 
the high-ionizati on [NeVI] 7.65 /im line, which can be 
bright in AGNs (jSturm et al. I 120021) and is impossible 
to separate from the PAH feature in low-resolution IRS 
spectra. Moreover, the 7.7 /im feat ure might be signif- 
icantly affected by AGN radiation (|Smith et al.l l2007bt 
IQ'Dowd et aT1l20M) . 

4.2. AGN dominated (clean sample) vs. starburst 
dominated sources 

In agreement with previous studi e s (e.g., 
Weedman et al.l 120051: iBuchanan et all 120061: iShi et al.1 
20061: iDeo et all 120091 : I Wu et al.l 120091 ) . we observe a 



large variety of infrared spectral shapes among Seyfert 
galaxies. However, those spectra dominated by dust re- 
processed emission generated by black hole accretion (see 
examples in the left column of Fig. [1]), rather than by star 
formation, almost invariably peak (when plotted in vF v 
units) at ^15-20 /im, in good agreement with models 
of dust y tori heated by a central source o f UV radiation 
(e.g.. iDullemond fc van Bemmell 120051: iHonig et al.l 
[20061: iNenkova et al.l l2008t lAlonso-Herrero et alj 120111 )7 
Furthermore, the AGN components of those IRS spectra 
with inferred significant starburst contamination (see 
the middle column of Fig. [T]) prove to be similar to 
the spectra of "pure" AGNs. In particular, most of 
the former also peak at 15-20 /im. However, since our 
procedure of estimating the star formation contribu- 
tion based on the strength of PAH features becomes 
progressively less reliable with increasing wavelength, 
there is much uncertainty in the deduced AGN spectral 
contributions at A > 20 /im. All these findings are 
similar to the results of previous attempts to decompose 
Spitzer spectra of quasars and Seyfer t galaxies into AGN 
and starburst comp onents (e.g., iNetzer et al.1 120071 : 
IMullanev etld1l20Tir) . 

The IRS spectra of six objects, NGC 1142, ESO 005- 
G004, IGR J07563-4137, NGC 4945, IGR J14561-3738, 
and MCG+04-48-002 (see the right column of panels in 
Fig.[TJ closely resemble the starburst template. We found 
it practically impossible to distinguish AGN and host 
galaxy components in these starburst dominated sources 
and therefore excluded them (see Table [l} from most 
of the subsequent analysis. Interestingly, most of our 
starburst dominated objects are strongly X-ray absorbed 
AGNs (Ah ~ 10 24 cur 2 ). One may speculate that i) 
large supplies of cold gas and dust associated with star- 
burst activity in a galactic nucleus facilitate the forma- 
tion of a dense central obscuring torus, and/or ii) part of 
the X-ray absorption is caused by cold gas tracing star 
formation in the galaxy and located outside a parsec- 
scale AGN torus. 

The remaining 61 objects (with NGC 1068 excluded 
for being a Compton-thick source) compose a "clean" 
sample for our subsequent analysis. As concerns the 
LLAGN NGC 4395, since its (low signal to noise) low- 
resolution IRS spectrum leaves doubts as to the presence 
of a significant starburst contribution (we estimate it at 
~ 14% at 15 /im), we have also analyzed available high- 
resolution IRS data, fol l owing the methods described aby 
iGoulding fc Alexander! (|2009j ). While the derivation of 



an accurate continuum shape, using only high-resolution 
IRS, is complicated by the tip-tilt effects of the individ- 
ual echelle orders, the high-resolution MIR spectrum of 
NGC 4395 is characterized almost entirely by an AGN- 
produced broken power law with little or no evidence 
(EW <C 0.1 /im) for superposed PAH features. Hence, 
NGC 4395 is clearly AGN dominated and thus should be 
part of our clean sample. Nevertheless, as noted before, 
we still distinguish this "dwarf Seyfert" from the rest of 
the sample during our HX-MIR cross-correlation analy- 
sis because it might represent a physically different class 
of AGNs. 

4.3. 15 fj,m flux and luminosity 

We have just seen (Fig. Q] and ^4.21) that after sub- 
traction of the star formation contribution, AGN MIR 
continua have an approximately constant shape. This 
suggests that it should be possible to estimate bolomet- 
ric luminosities of AGN obscuring tori using their flux 
densities measured at a single MIR wavelength. We have 
chosen to use for this purpose the rest- frame A = 15 /im. 
Specifically, /„ (15 /im) was determined by averaging a 
given spectrum over the wavelength range 14.7-15.2 /im. 
There are several reasons behind this choice. First, 
A = 15 /im is approximately where AGN torus emis- 
sion peaks. Second, since IRS SL2 data are not available 
for some of our sources, we can only use wavelengths 
A > 8 /im for the whole sample. Third, wavelengths 
A > 20 /im are disfavored because cool dust emission 
associated with star formation becomes more important 
with increasing wavelength and frequently dominates far- 
infrared and even mid-infared spectra of Seyferts. Fi- 
nally, there are no strong emission lines or absorption 
features within ~ 0.5 /im on either side of 15 /im. 

The last three columns of Table Q] present the total 
measured flux densities (/„) and corresponding luminosi- 
ties (yhv) at 15 /im as well as estimated fractions of AGN 
emission in the total flux at 15 /im, -Fagn ■ Statistical 
uncertainties for the infrared fluxes and luminosities are 
negligibly small. For the six starburst dominated ob- 
jects, we assume that -Fagn < 50%. In the subsequent 
analysis, the starburst subtracted MIR flux of AGNs is 
defined as 

/i5/«n = F AG ^vf v (15 /im). (1) 
The corresponding AGN luminosity is defined as 

£i5A«n = F AW vL v (15 /im). (2) 

4.4. Uncertainty in AGN MIR flux estimates 

The main potential source of systematic uncertainty 
in our estimation of AGN fluxes at 15 /i m is the use 
of a fi xed starburst spectral template from lBrandl et aLl 
(2006). In reality, the spectral properties of infrared 
emission from dust heated by starbursts may vary 
from one galaxy to another, and some authors (e.g., 
IMullanev et al.l 1201 ID have attempted to take this di- 
versity into account in separating AGN and host galaxy 
components for Seyfert galaxies. 

A crude estimate of the systematic uncertainty asso- 
ciated with our use of a fixed spectral template was al- 
ready made in tj4. II using the difference in the normaliza- 
tions of starburst spectral components determined using 
the 6.2 and 11.3 /im PAH features. Namely, starburst 
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component amplitudes could be estimated by our fitting 
procedure to within ~ 20%. This implies that for AGN 
dominated sources, i.e., those objects with Fagn ^ 50%, 
the AGN fluxes at 15 /jm are also estimated to within 
~ 20%, i.e., to better than 0.1 dex in log space. 

To better understand the uncertainties associated with 
our estimates of /ls^m, we performed an alternative 
spectral analysis based on a set of starburst spectral 
templates for a subsample of our objects. Specifi- 
cally, we selected a representative subset (16 objects) 
of AGN-dominatcd, mixed, and starburst-dominated 
sources. This comparison sample is equally divided be- 
tween sources with IRS-staring and mapping data. In 
H4.ll we assumed that the SL — LL discontinuity ob- 
served in the MIR spectra of some of these AGNs 
arises due to extended host galaxy emission which pro- 
vides an additional contribution to the larger aper- 
ture LL spectrum. However, it has also been sug- 
gested in previous studies that this discontinuity be- 
tween the spectral orders derives from an IRS detec- 
tor effect; for the purposes of comparison, we imposed 
this assumption for our sub-sample. Furthermore, for 
those sources with IRS mapping data, we combined 
the rogue-pixel cleaned BCD images and extracted nu- 
clear spect ra using the 3D sp ectral reduction program 
CUBISM (jSmith et all l2007al) . whic h is used widely 
in recent MIR AGN l i teratu r e (e.g.. iDale et al.l 12009b 
Gould mg fc Alexander! |2009t IDiamond-Stanic fc Riekd 
2010HPetric et al.ll20lfl(Alonso-Herrero et al.ll2012h . We 
used DecompIR (jMullanev et al.ll2010f ) to deconvolve the 
MIR spectra for our comparison sample using a set 
of empirical and the oretical starburst templates (e.g., 
iGoulding et al.l [201 ID and assumed an absorbed broken 
power law to model the AGN component. The derived 
AGN fractions and fluxes prove to be entirely consis- 
tent with the values established by our baseline analysis 
f q3.il N4- 1|) for those sources with Fagn > 0.5 (the aver- 
age scatter in the derived /i5 Mm is ±8%). However, for 
those sources with -Fagn < 0.5 (i.e., the spectra appear 
starburst dominated), the spectral fits become strongly 
dependent on the imposed starburst templates, and the 
scatter in the measured AGN flux increases to a factor 
of - 2. 

As an additional check of our IRS spectral mea- 
surements, we can use the IRAC imaging data avail- 
able for all of our AGNs. For the vast majority of 
the sources, the absolute values of flux densities mea- 
sured by IRAC and IRS in the overlapping spectral re- 
gion below 8 (im are in good mutual agreement (see 
Fig. [TJ, especially when the smallest (2.4") IRAC aper- 
ture is used (recall that the IRS SL slits have similar 
widths, 3.6-3.7"). However, there are a few sources 
for which there is a significant discrepancy between the 
spectroscopic and photometric fluxes at A < 8 /im. 
Since the corresponding IRAC and IRS observations 
were separated by several years in time, these flux dif- 
ferences probably indicate significant intrinsic variabil- 
ity of A GN torus emission, esp e cially at shorter wave- 
lengths dSuganuma et all 120061: IKishimoto etaLl 120091: 
iTristram et all 120091: IKozlowski et al.l 12010ft . Also, as 
shown in Appendix the detection of significant ex- 
tended emission at 8 fim by IRAC in many sources is fully 
consistent with our conclusions about the host-galaxy 
contamination of IRS spectra. 



Finally, we compared our results for 16 AGNs 
with higher angu lar resolution observations from 
iGandhi et all (|2009ft . We find excellent agreement (see 
Fig. [5] and discussion in £|5.2[) between the fluxes derived 
from the very different observational data sets, which 
suggests that host galaxy contamination of our measured 
values of AGN MIR fluxes is minimal. 

We conclude that the combined statistical and system- 
atic uncertainty in our AGN 15 /im-flux estimates for the 
clean sample (i.e., AGN dominated sources) is probably 
less than 0.1 dex. This uncertainty proves to be small in 
comparison with the intrinsic scatter in the HX-MIR flux 
and luminosity correlations (see below) . We further 
discuss the potential influence of starburst contamination 
on our derived correlations in £15.11 

5. HX-MIR LUMINOSITY RELATION 

Before comparing the luminosities of AGN struc- 
tural components (accretion disk, corona, and obscur- 
ing torus), which will be the subject of the next sec- 
tion (fcj6]), we first perform a cross-correlation analysis of 
AGN luminosities measured in the HX and MIR bands 
by INTEGRAL and Spitzer, respectively. The results of 
this analysis may also be interesting in its own right for 
any studies addressing links between X-ray and infrared 
emission in AGNs. 

Figure [5] shows the scatter plot of £i5 Mm vs. £<hx, 
where Lhx is the luminosity in the 17-60 keV energy 
band and Li5 pm was defined in eq. ©. In computing 
luminosities from fluxes, we neglected uncertainties as- 
sociated with source distances. 

Considering the clean AGN sample without NGC 4395 
and fitting i-is^m as a function of Z-hx (computing the 
linear regression in log- log space) , we find a strong, non- 
linear correlation between MIR and HX luminosities (see 
Tabled]): 

(1.7±0.2)L£™, (3) 

where the luminosities are measured in units of 
10 43 erg s _1 . The rms scatter of L 15flm values around 
the mean trend is 0.34 dex. 

As can be seen in Fig. (2J NGC 4395 is a clear outlier 
from the luminosity correlation, with its MIR luminos- 
ity being almost 2 orders of magnitude below the Lux~ 
-F15 trend described by eq. ([3]). If we consider this 
LLAGN together with the rest of the clean sample, the 
slope of the correlation increases from 0.74 to 0.85 (see 
Table O, although formally the change is not significant. 

Fig. [2] also shows the six starburst dominated sources, 
assuming that their AGN fractions Fagn < 50%. Sur- 
prisingly, the upper limits to the 15 /xm fluxes of AGN 
components for all these objects lie below the best-fitting 
relation for AGN dominated sources (eq. [3]). At least 
in some cases, this behavior is likely caused by attenu- 
ation of the intrinsic MIR emission from the nucleus in 
the obscuring torus a nd in the surrounding galaxy (see 
IGoulding et al.l [20121 1. In particular, the IRS spectrum 
of NGC 4945 (Fig. Q} exhibits a very deep silicate ab- 
sorption trough at 10 /im, which, assuming the standard 
composition of interstellar dust (Drain j |200l and the 
simplest scenario of an infrared source surrounded by a 
shell of dust, suggests that the neighboring continuum 
emission at — 15 /jm should be attenuated by a factor of 
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Table 2 

Results of HX-MIR cross-correlation analysis for the clean sample of AGNs and its subsamples: 



J15Min,43 



Sample 


Number 


a 


b 


rms, 




Spearman 




Pearson 




of objects 






dex 


P 


^"null 


r 


-Pnull 


Without NGC 4395 


60 


1.7 ±0.2 


0.74 ±0.06 


0.34 


0.85 


7 X 10- 18 


0.85 


8 x 10~ 18 


All 


61 


1.4 ±0.2 


0.85 ± 0.06 


0.39 


0.86 


9 x 10~ 19 


0.88 


3 x 10- 20 


Fagn > 0.9 


46 


1.7 ±0.3 


0.72 ±0.07 


0.35 


0.82 


5 x 10" 12 


0.84 


5 x lO" 13 


Without NGC 4395, total MIR fluxes 


60 


2.0 ±0.2 


0.68 ±0.06 


0.35 


0.83 


1.2 x 10~ 16 


0.82 


8 x 10~ 16 


z < 0.02 (without NGC 4395) 


35 


1.8 ±0.2 


0.93 ±0.10 


0.29 


0.84 


3 x lO" 10 


0.85 


1.3 x lO -10 


z > 0.02 


25 


1.7 ±0.8 


0.69 ±0.16 


0.38 


0.67 


2 x 10~ 4 


0.66 


3 x 10~ 4 


Syls and NLSyls 


33 


2.1 ±0.4 


0.69 ±0.09 


0.37 


0.81 


1.2 x 10~ 8 


0.80 


2 x 10" 8 


Sy2s (without NGC 4395) 


27 


1.7 ±0.2 


0.85 ±0.09 


0.30 


0.85 


3 x 10~ 8 


0.88 


1.7 x 10~ 9 




Figure 2. Luminosity scatter plot of L15 llm vs. Lhx- Filled 
circles represent AGNs from the clean sample (excluding six star- 
burst dominated AGN). The black solid line shows the best-ftting 
power law (L15 um as a function of Z/hx) f° r these objects exclud- 
ing NGC 4395 (eq. [3]), while the two black dashed lines show this 
dependence multiplied and divided by 2.19, the rms scatter around 
the mean trend. The black long-dashed line shows the best-fitting 
relation for the total clean sample including NGC 4395. The cor- 
relation parameters are listed in Table [2] Empty squares denote 
the 6 starburst dominated sources, for which an AGN fraction 
Fagn < 50% at 15 /an is assumed. Also shown is the Compton- 
thick Seyfert 2 galaxy NGC 1068. The magenta dotted line shows 
the resu lt of fitting Lhx as a function of L15 Mm by a power law, 
cq. (HI), for the clean sample excluding NGC 4395. The ma- 
genta dash-dotted line s how s the same dependence corrected for 
the Malmquist bias, eq. <B2t . 

~ 3-5. In reality, depending on the actual distribution 
of dust in the nucleus and body of the galaxy, the AGN 
MIR emission can be absorbed even more strongly than 
suggested by the depth of the 10 itm trough. 

Finally, Fig. [2] shows the Compton-thick Seyfert 2 
galaxy NGC 1068, which is a clear outlier from the cor- 
relation between Lhx and £15 Mm described by eq. ([3]). 
This result is expected because intrinsic hard X-ray emis- 
sion is strongly absorbed in this object, and it is only 
the infrared signal that reveals the true power of this 
AGN. In fact, the discrepancy between the general trend 
and the position of NGC 1068 on the Lnx~Li 5flm dia- 



Figure 3. Li5 mii1 /Lhx ratio as a function of Lhx for the clean 
sample of AGNs. The uncertainties associated with the luminosi- 
ties are less than ~ 0.1 dex. Averages over 0.5 dex-wide bins in 
Ljjx are a l so shown, with the vertical error bars illustrating the 
rms scatter of individual measurements within bins. 

gram suggests that its true X-ray luminosity is two or- 
ders of magnitude higher than measured by INTEGRAL, 
i.e., L H x, unabsorbcd ~ 10 44 erg s" 1 . This estimate is 
consistent with values reported in the literature (e.g., 
iMatt et alj|2000l) . 

The mean trend described by eq. ([3]) suggests that 
the MIR/HX luminosity ratio decreases with increasing 
-Lhx- This can be better seen in Fig. [3l which shows 
the Li5 Min /^HX ratio as a function of I/hx- Grouping 
our clean sample into 0.5 dex-wide bins in Lhx shows 
that the Lux ratio decreases from ~ 1-3 at 

Lux ~ 10 42 -10 43 erg s" 1 to - 0.3-1 at L H x ~ 10 44 - 
10 45 erg s-\ although the "dwarf Seyfert" NGC 4395 is 
a clear outlier from this trend. 

5.1. Robustness of the correlation 

The derived iHX _ -^i5^m relation, eq. ([3]), makes it pos- 
sible to predict the HX luminosity for a given MIR lu- 
minosity. The combination of three facts, i) that our 
AGN sample is hard X-ray selected, ii) that this sample 
is not limited by sensitivity in the MIR band, and iii) 
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that the correlation has been derived by fitting £15 pm as 
a function of Lhx, ensures that this relation reproduces 
the intrinsic correlation between Lhx and £15 Mm for the 
local AGN population without any bias. 

To further test the robustness of the derived trend, we 
repeated our Lhx^is cross-correlation analysis for 
various subsamples of AGNs. First, one may ask whether 
our procedure of separating AGN and starburst spectral 
components significantly affects the results. To address 
this issue, we computed the correlation for 46 strongly 
AGN dominated sources - those objects in which the 
AGN component accounts for at least 90% of the total 
emission at 15 fim (i.e., Lagn > 0.9). The result (see Ta- 
ble [2]) is very close to the correlation found for the total 
clean sample from which NGC 4395 is excluded (eq. [3]). 
In addition, we repeated the analysis for the clean sample 
(without NGC 4395) using total measured 15 /im fluxes 
instead of AGN fluxes (i.e., setting Lagn = 1)- The am- 
plitude of the correlation increased by ~ 20%, obviously 
due to the unsubtracted contribution of starburst emis- 
sion, but the slope changed by less than la from 0.74 to 
0.68. These tests demonstrate that the correlation be- 
tween Lhx and Li5 Mm is not significantly affected by the 
details of our spectral analysis of IRS data. 

We next repeated the analysis separately for nearby 
(z < 0.02, 35 objects, excluding NGC 4395) and distant 
(z > 0.02, 25 objects) sources from the clean sample. The 
correlation, in particular the slope of 0.69 ±0.16, derived 
for the distant subsample (see Tabled is fully consistent 
with the correlation found for the total sample (eq. [3]). 
Since the z > 0.02 set, owing to the INTEGRAL detec- 
tion limit, is represented by luminous AGNs only, with 
Lhx ~ 10 43 -10 45 erg s _1 , this result also implies that 
the slope of the high-luminosity part of the HR-MIR 
correlation is not significantly different from the trend 
found over a broader range of luminosities. However, the 
slope, 0.93 ± 0.10, determined for the nearby (z < 0.02) 
subsample, mainly consisting of lower luminosity AGNs 
with Lhx ~ 10 42 -10 44 erg s _1 , is somewhat different 
from the general trend, but this difference is less than 2a 
significant. 

Finally, we repeated our analysis for different types of 
AGNs, namely Seyfert Is (Syl, including the interme- 
diate types 1.2 and 1.5) and Seyfert 2s (Sy2, including 
the intermediate types 1.8 and 1.9). The derived rela- 
tions (Table [5]) are consistent with each other and with 
eq. ©. As can be seen from Fig. |U Syls and Sy2s do 
not distinguish themselves on the -C/HX-^iS/jm diagram, 
nor do narrow-line Seyfert 1 galaxies occupy a distinct 
region of this diagram. Finally, there is no significant 
dependence of the Lis^m/Lnx ratio on the X-ray ab- 
sorption column density except for the clear separation 
of the extremely Compton-thick source NGC 1068 from 
the rest of the sample. 

We conclude that the correlation between HX and MIR 
luminosities described by eq. ([3]) is robust, although there 
is a weak indication that the slope of the correlation is not 
constant and decreases with increasing AGN luminosity. 

For some applications, one may also be interested in 
knowing the distribution of Lhx for a given Li5^ m , 
rather than Li5 Mm as a function of Lhx- We have there- 
fore also computed (Appendix [Bj the inverse linear re- 
gression, i.e., log Lhx as a function of logLi 5Aim , for our 
clean sample excluding NGC 4395. As can be seen in 
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Figure 4. Scatter plot of Li5 Mm vs. Lhx for the clean sample 
of AGNs excluding NGC 4395. Blue filled circles and red squares 
denote Seyfert Is and Seyfert 2s, respectively. Also indicated are 
narrow-line Seyfert Is (empty circles) and the X-ray bright opti- 
cally normal galaxy IGR J13091+1137 (empty square). The solid 
black, dotted blue, and dashed red lines show the best-fitting power 
laws for the clean sample (without NGC 4395), Syls, and Sy2s, re- 
spectively (Table [2}. 

Fig. [2J this relation is different from the dependence of 
Li5fj,m on Lhx- 

5.2. Comparison with previous work 

iMullanev et al.l (|2011t ) have studied infrared proper- 
ties of nearby AGNs detected in the 14-195 keV energy 
band by Swift/BAT. This sample, although not statis- 
tically complete, is similar to our INTEGRAL sample 
in that it is hard X-ray selected. Using a subsample of 
AGNs having bot h Spitzer/IRS spectros copic and IRAS 
photometric data, IMullanev et al.l ([201 ID developed and 
tested a procedure, based on a set of starburst spectral 
templates, that make it possible to separate AGN and 
starburst contributions to the infrared flux using IRAS 
four-band photometry only. They then applied this pro- 
cedure to a sample of 44 BAT AGNs and found that 
Li2 Mm ,43 = (2.4 ± 0.4)L; 4 74 ± 5 1 1 f eV 4 3. This result is in 
excellent agreement wit h our eq . (|3|) . 

On the other hand, 1 Gandhi et al.l (|2009|) have re- 
ported a near proportionality between 2-10 keV (Lx) 
and 12 /im luminosities for Seyfert galaxies using high 
angular resolution infrared observations: their best esti- 
mate is Li2 A , m oc L x 11±0 07 . This result seems to contra- 
dict our conclusion that the Li5 Mm /LHx ratio decreases 
with increasing luminosity. 

A number of factors might contribute to this discrep- 
ancy, but the most important one appears to be the dif- 
ference in sample luminosities. A difference in galaxy 
weighting makes at most a minor difference. In their 
preferred regression procedure. lGandhi et alj (|2009f) took 
into account individual uncertainties in X-ray and in- 
frared luminosities. However, the X-ray uncertainties 
were estimated in a rather arbitrary way, taking into ac- 
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count long-term variability for some sources but not for 
others. This led to significantly different weights ascribed 
to different sources in fitting. In our view, for the prob- 
lem at hand, it is preferable to use a standard linear re- 
gression procedure in log-log space giving equal weights 
to all the sources in a sample. In fact, iGandhi et al.1 
( 2009) did perform such an analysis and obtained a some- 
what flatter dependence Li 2lim °c L^ 02±0 ' 07 , which is, 
however, still significantly steeper than the £15 Mm oc 
L a74±o.06 relation found here _ 

Another potentially importa nt factor is th e use o f 2- 
10 keV X-ray luminosities by IGandhi et al.l (|2009ft vs. 
our use of HX luminosities. Due to our hard X-ray selec- 
tion and the resulting insensitivity to absorption effects, 
we are able to determine HX luminosities directly from 
17-60 ke V fluxes measu r ed by INTEGRAL. In compar- 
ison, the IGandhi et al.1 (120091 ) sample contains a large 
number of significantly absorbed (N-r > 10 23 cm~ 2 ) 
sources whose intrinsic 2-10 keV luminosities were esti- 
mated through model-dependent analysis of X-ray spec- 
tra or, in some cases, even using [OIII] optical line 
fluxes. We have compared the HX and X-ray luminosi- 
ties for 16 Seyfert galaxies (excluding the Compton thick 
Seyfert NGC 1068) which a r e pres ent in both the IN- 
TEGRAL and IGandhi et all (|2009| 1 samples. The data, 
spanning the Lhx range from 9 x 10 41 erg s^ 1 (Cen A) 
to 1.4 x 10 44 erg s _1 (Mrk 509), prove to be consistent 
with Lux being proportional to Lx- The small scatter 
(0.23 dex) associated with this correlation can be fully 
attribute d to the uncertainties in the Lx values as esti- 
mated by IGandhi et al.l (|2009l ) , whereas the mean ratio 
Lhx/Lx ~ 1-5 may be considered typical for AGNs and 
is only slightly larger than the ratio (1.2) correspond- 
ing to a fiducial AGN spectrum used in our analysis be- 
low ( £16.11) . Furtheremore, sys tematic studies base d on 
AGNs from the INTEGRAL fie Rosa et al1 l2Tjrl and 
Swift/BAT (|Winter et al.ll2009h hard X-ray surveys have 
not revealed a significant dependence of the HX/X-ray 
flux r atio on luminos i ty. W e thus conclude that the use of 
L x bv lGandhi et~a l. (2009) vs. our use of Lux is unlikely 
to lead to a significant difference between the results of 
the corresponding X-ray-infrared cross-correlation anal- 
yses. 

Differences in the infrared data analysis might also 
play a role. Our study is based on Spitzer/IRS spec- 
troscopy and uses a template base d separation o f AGN 
and starburst spectral components. IGandhi et al.l ([2009D 
use narrow-filter photometry near rest- frame 12 fim with 
the VISIR instrument on VLT, which provides a signifi- 
cantly better angular resolution compared to Spitzer/IRS 
and hence presumably minimizes host galaxy contami- 
nation. We have again used the overlapping sample of 
16 Seyferts to check if there are any systematic differ- 
ences associated with these two approaches. Specifically, 
we used our IRS spectra and der ived fluxes within th e 
VISIR filters that were used by IGandhi et all ((2009), 
which differ from object to object. As shown in Fig. [Si 
although the total fluxes measured by IRS spectra are 
somewhat higher than those measured by VISIR, our 
standard correction for the host galaxy contamination 
brings both data sets to nearly perfect agreement. A sig- 
nificant difference (a factor of ~ 2) is only observed for 
NGC 3081, but it seems natural to expect some devia- 
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Figure 5. Co mparison of fluxes measured at A ss 12 fim 
by VLT/ VISIR UGandhi et al.ll2009h and Spitzer/IRS. The black 
empty squares are total IRS fluxes, and the black dots are 
star burst-subtracted fluxes, to which systematic uncertainties of 
20% are ascribed. The red cross with the arrow is the total flux for 
ESO 005-G004, for which IRS cannot confidently resolve the AGN 
component. The 1:1 dependence is shown with the solid line. 
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Figure 6. Distribution of AGN s from our clean sample (solid 
line) and the IGandhi et al.l (2009) one (dashed line) over the HX 
luminosity. For the latter sample we assumed that Lhx = 1.5Lx- 



tions in a sample of 16 objects, given that Spitzer/IRS 
and VLT/ VISIR observations were not simultaneous (for 
the same reason, the 8 /im fluxes measured by IRAC and 
IRS differ significantly for some of our AGNs, see §4.4|) . 
Furthermore, there is no trend with either distance or 
luminosity, although the comparison sample spans dis- 
tances from 3.6 Mpc (Cen A) to 174 Mpc (ESO 209- 
G012). We conclude that the differences in the infrared 
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data analysis do n ot significantly bias our results with 
respect to those of lGandhi et al.l (12009ft. 

Perhaps most importantly. IGandhi etaLl (|2009l ) used 
a heterogeneous sample of 42 Seyfert galaxies, whereas 
we use a statistically complete and somew hat larger sam- 
ple. As a result, the IGandhi et a l. (2009) sample is sig- 
nificantly shifted to lower luminosities relative to ours 
(Fig. [6]): e.g., the corresponding fractions of AGNs with 
£hx > 10 44 erg s" 1 are 12% and 38%. As was di- 
cusssed in M5 - 11 dividing our sample into two subsets, 
z < 0.02 and z > 0.02, represented by relatively low 
and high luminosity AGNs, respectively (~ 10 42 -10 44 vs. 
~ 10 43 -10 45 erg s _1 ), tentatively suggests that the slope 
of the X-ray-infrared correlation changes from 0.93±0.10 
to 0.69 ± 0-1 6 as t he AGN luminosity increases. The 
IGandhi et al.l (|2009l ) sample effectively probes the lumi- 
nosity range Lhx ^ 10 44 erg s _1 , similar to our z < 0.02 
subsample, and the slopes inferred for these two data 
sets are in satisfactory agre ement with on e anoth er. This 
suggests that the results of lGandhi et al.l (|2009D and the 
present work are actually consistent with each other. 

We conclude that further studies using larger, well- 
defined samples of AGNs are required to clarify if the 
slope of the X-ray-infrared correlation depends on lumi- 
nosity, as tentatively suggested by the existing data. 

6. TORUS VS. DISK AND CORONA 

The unified model posits that a torus of molecular gas 
and dust subtending a solid angle Sl-porus (if viewed from 
the SMBH) intercepts optical, UV, and soft X-ray radi- 
ation from the central accretion flow and converts it into 
thermal infrared emission. Therefore, assuming that the 
central source of radiation is isotropic, the total luminos- 
ity of the torus is expected to be 



L 



n 



Torus 



Torus 
47T 



^Disk, 



(4) 



where LDisk is the luminosity of the accretion disk, pre- 
sumably emitted between A ~ 1 /im (NIR) and E ~ 
2 keV (soft X-rays). These boundaries usually separate 
the MIR, BBB, and HX components (see TO) in the SEDs 
of typ e 1 AGNs (see, e.g-. lElvis et al.lll994l:1Sazonov et al.l 
[200l . Physically, the 1 /im boundary marks the on- 
set of thermal emission from hot dust at sublimation 
temperature (~ 1500 K), whereas accretion disk emis- 
sion is expected to peak in the near - or far-UV bands in 
quasars and Seyfert gala xies (e.g., IShakura fc Sunvaevl 
H973tlHubenv et al.ll200l . Therefore, the chosen energy 
boundaries for Loisk ensure that virtually all of the ac- 
cretion disk luminosity is accounted for. 

Depending on the column density through the torus, it 
can also reprocess a fraction < Otoi-us/ (47r) of the higher 
energy (2-10 keV) luminosity emitted by a hot corona 
of the accretion disk. We have neglected this contribu- 
tion in eq. dU, first because the torus's optical depth 
may be significantly smaller than unity for 5-10 keV X- 
rays, in contrast to the softer emission from the accre- 
tion disk, and also because we expect the X-ray (below 
10 keV) luminosity of the corona to be much lower than 
the bolometric luminosity of the accretion disk. This 
last assumption will be verified below upon completion 
of our cross-correlation analysis. Finally, it is assumed 
that none of the hard X-ray emission (above 10 keV) 



is absorbed within the AGN, which is a reasonable as- 
sumption expect for very Compton-thick objects such as 
NGC 1068. 

As demonstrated below, Spitzer and INTEGRAL data 
together make it possible to estimate the luminosity 
(^Torus) and the solid angle (f^Torus) of the torus as well 
as the luminosity of the corona (at energies 2-300 keV), 
^Corona- We can therefore use eq. (|4]) to study the rela- 
tionship between Lnisk and ^corona, i.e., between emis- 
sion properties of the accretion disk and corona. 

6.1. Bolometric corrections 

We proceed by determining coefficients for conver- 
sion of the measured quantities Lhx and Li5 Mm to the 
AGN intrinsic quantities ^corona and Ltoi-us, respec- 
tively. Hard X-ray spectral shapes do not vary much 
from one Seyfert galaxy to another, apart from the pho- 
toabsorption rollover in type 2 objects below 10 keV. 
Typically, absorption corrected AGN spectra can be de- 
scribed above 2 keV as a power law wit h a photon in- 
dex r ~ 1.7 (e.g.. IReeves fc Turn er 2000) and a rollover 
above ~ 100 keV (e.g.. IMolina et all 120091) . We adopt 
that 

^HX « 0.3-Lcorona- (5) 

This relation corresponds to a power-law spectrum with 
T = 1.7 and an exponential cutoff at E{ — 200 keV 
and is consistent with an average 3-300 keV spectrum 
of - 100 local AGNs detec ted during INTEGRAL an d 
RXTE surveys of the sky (jSazonov et all [20081 [20Toh . 
Assuming that the values of the power-law index and 
cutoff energy vary from T — 1.5 to 1.9 and from Ef ~ 50 
to ~ 500 keV from one Seyfert g alaxy to another (as 
indicated by numer ous studies, e.g-. IZdziarski et aT1ll995l ; 
IMolina et al.M2009f ). we can estimate that the conversion 
described by eq. © can introduce a scatter in Lhx values 
for a given Lcorona of < 20%, i.e., < 0.1 dex. 

We next introduce a similar correction factor for the 
reprocessed emission from the torus: 



Li 



5 /-f.m 



0.5L To 



(6) 



To obtain the above coefficient, we compared the 
15 /xm luminosity with that integrated over the rest- 
frame 6-32 /im band, 1/6-32 ^mi f° r those Spitzer/IRS 
spectra (21 in total) that span this whole wavelength 
range (i.e., there are available data from the IRS second- 
order SL module) and do not suffer from significant 
contamination by MIR emission from dust associated 
with star formation. For the majority of these ob- 
jects, the ratio I/i5 M m/^6-32/jm is bounded in the narrow 
range of 0.65-0.85, only slightly depending on whether 
AGN emission lines (such as [OIV] 25.9 /im) are taken 
into account or not. We therefore adopted the relation 
Li5/j,m/Ije-32)j,m = 0.75 for AGN tori and additionally 
lowered this ratio by one third in eq. © to account for 
non-negligible (~ 50%) additional torus emission both 
shortward of 6 /xm and longward of 32 /im (see, e.g., 
INenkova et al.l [20081 ). While the resulting L i5 Mm /LTorus 
ratio (eq. [5]) is determined less strictly than the ratio 
Lax/ L corona above, the associated scatter in individual 
^15 fim/ -^Torus values around the mean trend given by 
eq. © is probably less than 20%, as suggested by the 
comparison of IRS spectra for "pure" objects, described 
above. 
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AGNs has a power-law dependence on luminosity: 



42 43 44 

Log L HX , erg s" 1 

Figure 7. Fraction of obscured (JVh > 10 22 cm" 2 ) AGNs in 
the local Universe as a function of hard X-ray luminosity, based 
on the INTEGRAL sample. The error bars represent the Poisson 
uncertainty associated with the number of objects in a given bin. 
The lowest and highest luminosity bins contain just one source each 
(NGC 4395 and IGR J09446-2636, respectively). The approximate 
description of the observed trend by eq. (1 1 3 I t is shown by the dashed 
line. 



6.2. Solid angle of the torus 

The next step in our analysis is to derive the torus 
solid angle ilxorus- To this end, we assume that, for a 
given hard X-ray luminosity, riibrus is determined by the 
relative number of obscured (type 2) AGNs of that lumi- 
nosity, i.e., 



n 



Torus 
47T 



ax) 



A^typc 2(£hx) 

AtotalC^HXj 



(7) 



We consider an AGN obscured if its X-ray absorption 
column density Ah > 10 22 cm~ 2 . In this connection, 
recall (see $2]) that the Ah values for our objects are 
not based on INTEGRAL hard X-ray measurements but 
have been determined through analysis of high-quality 
X-ray spectra obtained by various X-ray telescopes. 

Dividing our AGN sample into several bins in Lhx, 
Fig. [7] shows the fraction of absorbed AGNs as a func- 
tion of Lhx- A strong trend of decreasing At ypo 2/Aftotal 
ratio with increasing luminosity is evident. We can es- 
timate the statistical significance of this trend using a 
maximum-likelihood estimator, 



£ = -2^1nP 1>2 (L H x,i)- 



(8) 



Here, the summation is over our sample of AGNs and 
Pi ,2 is the probability for a given object with luminosity 
Lnx,i to be either obscured (P2) or unobscured (Pi). We 
restrict our consideration to the luminosity range 41.5 < 
loginx < 45, because there are only two objects in the 
sample which fall outside this range (one on either side, 
see Fig. [7]). Suppose now that the fraction of obscured 



Obviously, 



P 2 = p + a(logi HX -41.5). 



P = 1 - P 2 



(9) 



(10) 



Requiring that < Pi ,2 < 1 over the 41.5 < logi-nx < 
45 range yields constraints on the slope and intercept of 
the trend: 

0<p<l (11) 



and 



P 

3.5 



< a < 



1-P 
3.5 ' 



(12) 



With these definitions, the maximum of the likelihood 
function proves to be at p w 1 and a = 0.25. By in- 
tegrating C over (p, a) parameter space with the priors 
given by eqs. (fTTj) and (fl"2"]) . i.e., using a Bayesian ap- 
proach, we find that the probability that a < is 0.999. 
Hence, the declining luminosity trend of the obscured 
fr action is ascertained with more than 3er significance. 

iSazonov et al.l (|20100 have recently confirmed this lu- 
minosity dependence (see their Fig. 3) using a nearly 
doubled sample of AGNs detected during 7 years of 
INTEGRAL observations (compared to the 3.5-year 
all-sky survey used in the present work). Further- 
more, the existence of this trend has been reliably es- 
tablished in the past d ecade using vario us X-ray se- 
lected samples of AGNs (lUeda et a jjl2003t ISte ffen et alJ 



2003tlSazonov fc Revnivtsevll2q04l:lLa Franca et al.l200E 
Sazonov et al.l I2007t IHasingerl 120081: [Burlon et aljboilt 
se e in particular Fig 8 in IHasingerl 120081 and Fig. 13 
in iBurlon et al.l I2011D . Therefore, in accordance with 
eq. ([7]) we adopt the expression 



,1, logi H x<41.5 
»Torus = I 1 - 0.25(log L HX - 41.5), 
4?r 1 41.5 < logLnx < 45 
0.125, logLnx > 45. 



(13) 



We thus assume that the phenomenon of decreasing 
fraction of absorbed AGNs with increasing luminosity 
reflects an underlying trend of increasing opening angle 
of the obscuring torus. This is one of the crucial points 
in our analysis. According to eq. (|13[) . the slope of the 
^Torus (Lux) dependence is approximately equal to 0.25 
for Lux ranging between ~ 10 415 and 10 45 erg s _1 with 
the associated uncertainty being small, ~ 10% as de- 
termined from the dispersion of data points in Fig. [7] 
and from the Bayesian analysis described above. How- 
ever, the f^Torus (Lux) dependence holds true in a sta- 
tistical sense only, and there might be variations in the 
torus opening angle among AGNs of a given luminos- 
ity. Unfortunately, observations do not yet provide reli- 
able information on the distribution of f^Torus values for 
a given Lhx, and hence we cannot predict to what de- 
gree this scatter could affect our results below. Further- 
more, the exact parameters of the At ype 2 /A^ t a i (Lux) 
dependence adopted in eq. (fT3| should be applied only 
to the local (z ~ 0) population of AGNs, in particular 
because the fraction of obscured sources among high- 
luminosity AGNs ap pears to be larger in the distant 
(z > 1) Universe (e g. , iHickox et alll2007t lHasingerll200l 
iTreister et al|[2008h . 
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6.3. Disc vs. Corona 

We are now ready to estimate the accretion disk lumi- 
nosities for our AGNs using eq. (j4]): 



^Disk 



4?T 



Torus 



^Torus- 



(14) 



Specifically, we first determine Lcorona and Ltotos using 
eqs. (J5j> and ([6]), respectively and then use eq. (fl"4f to find 
i'Disk- The resulting scatter plot of LDisk vs. Lcorona f° r 
the clean sample excluding NGC 4395 is shown in Fig.[H 
Fitting the ^Disk vs. Lcorona data in Fig. [S] with a 
power law yields the correlation 



iDisk,44 = (1.59 ±0.16)L 



0.97±0.06 

Corona, 44' 



(15) 



where the luminosities are measured in units of 
10 44 erg s _1 . The rms scatter around the mean trend 
is 0.33 and 0.34 dex along the ^Disk and -^Corona coordi- 
nates, respectively. 

The derived relation, eq. (TT51) , allows one to predict the 
disk luminosity for a given coronal luminosity. Hence, 
if the coronal luminosity -Lcorona of an AGN is known, 
one can expect its accretion disk luminosity to be equal 
within a factor of « 2 (at the la confidence level) to 
1.6L Co rona- The L Disk / L Co rona ratio does not signifi- 
cantly depend on luminosity in the effective range of 
^Corona from ~ 10 43 to 10 45 erg s- 1 . 

The 2-10 keV energy band contains ~ 25% of the total 
coronal luminosity. If all of this X-ray emission were con- 
verted in the torus into infrared radiation as efficiently 
as accretion disk emission, it would increase L^ orus by 
only ~ 15%. This justifies the approximation adopted in 
eq. flU). 

For some applications, one may also be interested in 
knowing the distribution of ^Corona f° r a given ^Disk, 
rather than ^Disk as a function of ^corona- We have 
therefore also computed (Appendix [B]) the inverse linear 
regression, i.e., logLcorona as a function of logLDisk, for 
our clean sample excluding NGC 4395. As can be seen 
in Fig. [H this relation is different from the dependence 

of ^Disk on ^Corona- 

6.4. Scatter around the mean trend 

Fig. [9] shows the distribution of residuals for the 
-kcorona-^Disk correlation. Although it is plotted 
in terms of Lrjorona deviations, the distribution of 
(HogLDisk = log I/Disk - (logi Disk (Lcorona)) residu- 
als is quite similar. The distribution can be well- 
described by a log-normal function, dN/d log Lcorona <x 
cxp[— (<51ogLcorona) 2 /2er 2 ], where a = 0.34 is the mea- 
sured rms scatter of the correlation. 

Although the observed scatter in the correlation is 
fairly small, the intrinsic correlation between the corona 
and disk luminosities in Seyfert galaxies is probably even 
tighter. Indeed, the 15 /im fluxes measured by Spitzer, 
from which the ^Disk values were derived, presumably 
represent reprocessed accretion disk emission averaged 
over a number of years. Interferometric observations 
of Seyfert galaxies have established that the size of the 
MIR-emitting d ust region is of the or der of several light 
years or more (|Tristram et all 120091) . Specifically, the 
characteristic size of the 12 fim source was found to scale 
approximately as the square root of the AGN luminos- 
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Figure 8. Inferred luminosity of the accretion disk vs. that of the 
hot corona for the clean sample of 60 AGNs (the LLAGN NGC 4395 
is excluded). Various AGN types are indicated by different symbols 
as in Fig. [4] Also shown is the Compton-thick Seyfert 2 galaxy 
NGC 1068, which was excluded from the analysis. The black solid 
line is t he b est-fitting power law dependence I/Disk (^Corona) given 
by eq. 1151 . while the two black dashed lines show this dependence 
multiplied and divided by 2.14, the rms scatter of the correlation. 
The magent a do tted line is the best-fitting dependence Z/Corona 
(^Disk) ( e 1- IB30 . and the magenta dash-dotted line is the same 
dependence corrected for the Malmquist bias (eq. IB4I ). 
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Figure 9. Binned di strib ution of the residuals of the -Lcorona - 
Lrjisk correlation, eq. 1151 1. The error bars represent the Poisson 
uncertainty associated with the number of objects in a given bin. 
The solid line shows the log-normal distribution corresponding to 
the measured rms scatter (a = 0.34). 
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itjQ and to vary from I ~ 0.5 pc for L\ 2 Mm ~ 10 42 erg s 1 
to I ~ 50 pc for L12 um ~ 1 46 erg s" 1 . For exam- 
ple, | Tristram et all (|2009T > and iTristram fc Schartmannl 
H201J found I « 2-3 pc for NGC 1068, NGC 1365, 
MCG 5-23-16, and NGC 4151 and ~ 10 pc for NGC 7469. 
Therefore, the LDisk values used in our cross-correlation 
analysis should represent accretion disk luminosities av- 
eraged over the ~ 2//c-long period immediately preced- 
ing the Spitzer observation of a given AGN, which is ex- 
pected to range from a few years for the least luminous 
sources to a few tens of years for the most luminous ones. 
As demonstrated in Appendix [C] variability of the HX 
coronal emission, detected by INTEGRAL on time scales 
shorter than the characteristic time scale of MIR varia- 
tions, is expected to induce a scatter of ~ 0.2-0.25 dex 
around the mean -E/Corona-^Disk trend. 

Additional contributions to the observed scatter 
around the mean £c rona _ E'Disk trend can be provided 
by systematic uncertainties associated with: i) measure- 
ment of Lux and Lis^m, each < 0.1 dex (see Table [T] 
and fl4.4p . ii) conversion from Lhx to Lcorona and from 
Lispm to i-ibms, also < 0.1 dex each ( §6.ip . and hi) de- 
termination of the mean torus solid angle (r^Torus) for a 
given AGN luminosity and consequently conversion from 
-^Torus tO iDisk (via eq. [H]), < 0.05 dex (g^U). Hence, 
each of the above effects can contribute of the order of, or 
less than, 0.1 dex to the observed scatter in the Lcorona~ 
Lnisk correlation. Adding these contributions in quadra- 
ture to that expected from varibility implies that the 
total induced scatter is ~ 0.25 dex. After subtraction 
of this contribution from the measured scatter of the 
icorona-^Disk correlation, with rms = 0.33 dex, there 
remains a scatter ~ 0.2-0.25 dex, i.e., a factor of 1.5-2. 

Therefore, the intrinsic correlation between accretion 
disk and coronal emission in Seyfert galaxies is fairly 
tight. Furthermore, neither the icorona-^Disk relation 
nor the Lnx~Li5 fim relation, from which it originates, 
shows a significant dependence on either optical AGN 
type or X-ray absorption column density, although there 
are exceptions, which are discussed below. Assuming 
that the torus is co-aligned with the accretion disk and 
is a quasi-isotropic MIR emitter, this suggests that the 
coronal hard X-ray emission is at most modestly (less 
than a factor of ~ 2) anisotropic. This conclusion 
also holds true if the obscuring torus is oriented quasi- 
randomly with respect to the accretion disk/corona axis 
because for given Li^^ m there is little scatter in Lhx- By 
the same argument, the hard X-ray luminosity of AGNs 
cannot be dominated by collimated emission from rela- 
tivistic jets. 

As to the origin of the remaining (unaccounted for) 
scatter, at least two effects are likely to contribute to 
it. First, our analysis was based on the assumption that 
the characteristic solid angle subtended by the obscuring 
torus, fiibrus, i s the same for all AGNs of a given lumi- 
nosity. In reality, it is po ssible that f^T orus varies from 
one object to another (e.g.,[Elitzur 2012j) . This would di- 
rectly affect our estimates of -t-Disk from Librus and intro- 
duce scatter in the resulting correlation between Lcorona 
and LDisk- Similarly, the amplitude of the Compton re- 

9 This is consistent with the simple argument based on consid- 
cring dust heating by a central source of UV emission IBarvainisI 
[1983). 



flection component, ignored in our analysis, may vary 
from one AGN to another, which would introduce addi- 
tional scatter. Therefore, the scatter in the intrinsic ratio 
of powers generated in the accretion disk and corona is 
likely even smaller than 0.2 dex. 

6.5. Comparison with typical quasars 

The results of this work pertain to nearby Seyfert 
galaxies, and it is interesting to put them into the broader 
context of the cosmic history of SMBH growth. To this 
end, we compare our findings wit h the propertie s of th e 
SED of the "average quasar" from iSazonov et all (|2004ft . 
This template essentially rests on the assumption that 
the cosmic X-ray background (CXB) represents the in- 
tegrated hard X-ray emissio n of all AGNs in the Uni- 
verse and on the argument of lSoltanl (jl982D that the cu- 
mulative bolometric luminosity of AGNs is determined 
by the mean radiative efficiency e with which the in- 
tegrated mass of local SMBHs has bee n accumulated 
over t he cosmic time. Adopting e = 0.1, Sa zonov et al.l 
(2004) found that « 12.5% of the bolometric luminos- 
ity (below 300 keV) of the average quasar, is emitted at 
energies above 2 keV. Attributing this emission to the 
corona of the accretion disk and making a small correc- 
tion for absorption in the 2-10 keV energy band (be- 
cause here we are interested in intrinsic rather than ob- 
served properties of accretion disks and coronae), we find 
that iDisk ~ 6Lcorona for the average quasar. If we in- 
stead assume that e = 0.06, which corresponds to the 
standard Shakura-Sunyaev disk around a Schwarzschild 
black hole, then (iDisk/^Corona) ~ 3.5. By construction, 
this ratio primarily characterizes quasars with £c orona ~ 
10 44,5 erg s" 1 , w hich produce the bulk of the CXB (e.g., 
lUeda et al.l[200l . 

We can now determine the corresponding average ra- 
tio for the local AGN population, using a completely 
different method. The mean trend given by eq. (|15[) 
and the associated scatter (0.34 dex) imply that the 
iDisk/icorona ratio varies between w 0.75 and « 3.4 (the 
la range) around the mean value of ~ 1.6. Averaging 
over the log-normal distribution of -t<Disk/-E'Corona yields 

(iDisk/icorona) = / (-^Disk/ -^Corona) G^Disk/ J dioisk ~ 

2.1, independently of Lcorona- This ratio character- 
izes the summed radiation of the local AGN population. 
Therefore, the (-E-Disk/icorona) ratio appears to be larger, 
but only by a factor of ~ 2, for typical quasars making 
up the CXB relative to typical AGNs in the local Uni- 
verse. This implies that the ratio of the disk and coronal 
luminosities is approximately constant in all actively ac- 
creting, radiatively efficient SMBHs. 

6.6. Possible effect of anisotropic accretion disk 
emission 

Our treatment so far has been based on the assumption 
that the accretion disk is an isotropic source, i.e., its ob- 
served luminosity is independent from the viewing angle. 
This implies that the primary (UV) and reprocessed (in- 
frared) luminosities of AGNs are related through eq. (j4|) . 
However, the actual angular distribution of radiation 
emergent from the accretion disk around a SMBH might 
be close to Lambert's law, in which case the luminosity 
per solid angle c?LDisk/ 'd£l oc cos#, where 9 is the viewing 
angle with respect to the axis of the disk. If, in addition, 
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Figure 10. As Fig. [8] but for an accretion disk emitting according 
to the cosine law and a dusty torus lyi ng in the plane of the disk. 
The best- fitting relation is given by eq fl7p. 

an obscuring torus is not randomly oriented but its equa- 
tor lies in the plane of the accretion disk, then there will 
be a different relation between Lxorus and Loisk: 



J Torus 



Torus 
47T 



L D 



isk* 



(16) 



Considering this possibility a feasible alternative to 
our baseline scenario, we repeated our calculations us- 
ing eq. (fi"6l) instead of eq. (|4]). The resulting correla- 
tion between the corona and disk luminosities is shown 
in Fig. [TU] and the corresponding best-fitting relation is 
given by: 



L D isk,44 = (3.4±0.3)L 



1.21±0.06 
Corona, 44' 



(17) 



with the rms scatter a — 0.33 in LDisk- 

The anisotropic scenario predicts significantly larger 
-^Disk/icorona ratios at high luminosities (a factor of ~ 3 
at Lcorona = 10 45 erg s _1 ) compared to the isotropic 
case. However, at these luminosities there is large un- 
certainty in the inferred Loisk values for the anisotropic 
case, which is not fully reflected in the formal uncertain- 
ties quoted in eq. (fTT)) . Specifically, the uncertainty as- 
sociated with our adopted dependence of the torus solid 
angle on luminosity, eq. (1131) . is not taken into account. 
As is clear from Fig. [71 this additional uncertainty is large 
at high luminosities, Lhx ~ 10 44 er S s_1 > which was rela- 
tively unimportant when we considered the isotropic sce- 
nario but is the dominant source of uncertainty for the 
anisotropic scenario due to the quadratic dependence on 
^Tbrus in eq. {16]). 

The above consideration ignores two potentially impor- 
tant effects. First, radiative transfer in the accretion disk 
may lead to a limb-darkening effect, making the emission 
even more collimated in the polar direction than cos ff 
(e.g., iSunvaev fe Titarchukl [l98l lLaor fe Netzen [l989h . 
On the other hand, strong gravity in the vicinity of a 



black hole tends to harden the spectrum and increase 
the luminosity for observers viewing the accretion disk 
at grazing angles, with this effect being e specially pro- 
nounced for rapidl y rotating black holes (jCunninghaml 

Taking these various factors into account, it is likely 
that the correlations obtained in our isotropic (Fig. [5J 
and anisotropic (Fig. [TO]) scenarios bracket the true 
relationship between the corona and disk emission in 
AGNs. Furthermore, orientation effects and object-to- 
object variations in the black hole spin can contribute to 
the scatter in the relation between Loisk and Lrjorona- 

7. DISCUSSION AND CONCLUSIONS 

The main result of this work is that the luminosities of 
the accretion disk and corona are nearly proportional for 
nearby AGNs: L D i s k oc L°^±° 06 . To reveal this correla- 
tion, we derived intrinsic (presumably emitted between 
A ~ 1 /im and E ~ 2 keV) accretion disk luminosities 
from observed torus luminosities, LTorus, using eq. (|14p . 
which assumes that radiation from the accretion disk is 
reprocessed in a dusty torus whose effective solid angle 
decreases with increasing luminosity from almost 47r at 
L HX ~ 10 42 erg s" 1 to - 4tt/8 at L HX ~ 10 45 erg s -1 , 
as suggested by the observed falling fraction of obscured 
AGNs (eqs. [7] and [IS])- This effect of decreasing ob- 
scuration fraction is responsible for the MIR/HX lumi- 
nosity ratio decreasing with increasing AGN luminosity: 



L 



cx L 



0.74±0.06 



15 (im ^ ^HX 

The observed £corona - ^Disk relation implies a mean 
ratio (iDisk/icorona) ~ 2 for nearby AGNs. For com- 
parison, for typical quasars producing the cosmic X-ray 
background, (L Disk /LcoTona) ~ 3.5-6 (see N6.5p . Hence, 
hard X-ray radiation from accretion disk coronae (with a 
possible contribution from jets) carries a significant and 
roughly constant fraction, ~ 15-35%, of the bolomet- 
ric luminosity of SMBHs accreting in radiatively efficient 
mode (with a possible exception of black holes accreting 
near the critical Eddington rate, see j|7J 51 below ). 

7.1. Intrinsic ratio of disk and corona luminosities 

The measured LDisk/Lrjorona ratio likely overestimates 
the ratio of the intrinsic UV and hard X-ray luminosi- 
ties produced by the accretion disk and corona, respec- 
tively, because roughly half of the luminosity emitted by 
the corona is intercepted by the disk (and a small addi- 
tional amount by the obscur ing torus) and only ~ 10- 
20% of this radiation (e.g., lHaardt fe Maraschil [l993h 
is reflected, while the rest is reprocessed into thermal, 
softer emission and thus contributes to the disk lumi- 
nosity. In contrast, it is well-known that the corona in 
Seyfert galaxies intercepts only a small fraction of the 
disk' s radiation, probably because the corona is patchy 
(e.g. JZdziarski et aIll997D . Indeed, the measured shapes 
of the hard X-ray spectra of Seyferts imply that the 
hot corona is characterized by an amplification factor 
A ~ 5-10 (jGilfanov et al.ll2000D . i.e., the luminosity of 
the Comptonized hard X-ray radiation emergent from 
the corona is several times the luminosity of the incident 
soft photons. This would imply that LDisk "C Lrjorona if 
the disk were embedded in a homogeneous corona. Since 
in reality LDisk ^ Lcoronai as suggested by the results of 
the present study, the presence of a strong big blue bump 
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in th e spectra of type 1 AG Ns, and theoretical arguments 
(e.g.. lZdziarski et ai"1ll997l) . the solid angle subtended by 
the corona with respect to the accretion disk must be 
small, i.e., to a first approximation the corona does not 
shield the disk from the observer. 

On the other hand, the Compton reflection hump is lo- 
cated approximately in our working INTEGRAL energy 
band (17-60 keV) and may significantly contribute to the 
measured hard X-ray flux (see examples of tentatively 
detected re flection components in INTEGRAL spectra 
of AGNs in lde Rosa et al.ll2012f ). so that we likely over- 
estimate the flux of HX radiation coming directly from 
the corona in our simplistic analysis. 

Considering these counteracting effects together, the 
intrinsic ratio of the disk and coronal luminosities 
could be a factor of < 2 smaller than the measured 
^Disk/^Corona ratio. This implies that in Seyfert galax- 
ies, approximately equal powers are generated in the ac- 
cretion disk and hot corona. 

7.2. AGN bolometric corrections. Current rate of 
SMBH growth 

Our results suggest that hard X-ray luminosity is a 
good proxy for bolometric AGN luminosity (Lboi = 
-^Disk + ^Corona), except for extremely Compton-thick 
sources like NGC 1068. Using eqs. (fT5)) and (jSj) , we can 
estimate the bolometric correction for the 17-60 keV en- 
ergy band and the associated la range (due to the scatter 
in the Lcorona-^Disk correlation): 

«9 (6-15) (18) 

in the range Lhx ~ 10 42 -10 44 5 erg s _1 . 

The cumulative hard X-ray (17-60 keV) luminosity 
density of nearby AGNs found by integrating their lu- 
minosity function measured by INTEGRAL at Lhx ?S 
10 41 erg s" 1 is (1. 2 ± 2) x 10 39 erg s" 1 Mpc"^ 
()Sazonov et al.l 120071 feOlOfl . Low- luminosity (L H x < 
10 41 erg s _1 ) AGNs may add up to ~ 50% to this volume 
emissivity, as follows from cross-correlating the cosmic 
X-ray backgr ound intensity with the local galaxy distri- 
bution (e.g.. iRevniytsev et alj|2008t iMivaii et "all 119941: 
ICarrera et al.l I1995D . Using the iboi/^HX ratio from 
eq. (|18p , the bolometric luminosity density of local AGNs 
is thus ~ (1 — 3) x 10 40 erg s _1 Mpc~ 3 . This implies 
that the integrated present-day growth rate of SMBHs is 
Msmbh [z = 0) ~ (2-5) X 10; 6 (0.1/e) M Q yr" 1 Mpc" 3 , 
where e is the average radiative efficiency of accretion. 
Comparing AfgMBH (z — 0) with the total mass density 
of local SMBHs, prmrh (z = 0) ~ (3-5) x 10 5 M a) Mpc" 3 
(|Yu k, Tremainel 120021: IMarconi et ail I2004D . the total 
SMBH mass is currently growing on a time scale ~ 
10 times the age of the Universe. 

This estimate of the SMBH growth rate does not fully 
account for the contribution of obscured accretion tak- 
ing place in Compton-thick (Ah 3> 10 24 cm~ 2 ) AGNs. 
Moreover, it is valid only for accretion that is occurring in 
a radiatively efficient mode. In reality, a substantial frac- 
tion of SMBH growth at the present epoch may be tak- 
ing place through a radiatively inefficient mode of accre- 
tion, dominated by mechanical rather than radiative en- 
ergy output (e.g.. lChurazov et "ah ; 2005: M erloni fc Heinzl 



2008]). Therefore, the total SMBH accretion rate may be 
higher. 

We can re-calculate the bolometric correction given by 
cq. (|18p to the standard X-ray band, 2-10 keV. Assum- 
ing, as in ^6.11 a power-law spectrum with T = 1.7 and an 
exponential cutoff with E f = 200 keV, the 2-10 keV/17- 
60 keV luminosity ratio Lx/^hx = 0.82 and therefore 

^wll (7-18). (19) 
Lx 

This formula predicts the bolometric luminosity of an 

AGN from its intrinsic (i.e., unabsorbed) luminosity in 

the 2-10 keV energy band in the range Lx ~ 10 42 - 
10 44.5 erg s -i 

Finally, we can estimate the bolometric correction for 
the MIR band (A = 15 fim), using eqs. (HHJ) and ©: 



J^15 fan 

This formula should be accurate to within a factor of 
- 2 for AGNs with L lbfira ~ 10 42 -10 44 5 erg s" 1 . We 
have not tried to take into accout the non-linear depen- 
dence of Li5 Mm on Lhx, because it should only be used 
to predict £15 Mm for a given Lhx but not Lhx for a 
given £15 Mm . A more reliable bolometric correction for 
the MIR band could be obtained by using a MIR se- 
lected sample of AGNs. Interestingl y, ea. (l20l) is in good 
agreem ent with an early estimate bvlSpinoglio fc Malkanl 
(1989), £boi/Li2 M m ~ 5, based on direct integration of 
the IR-to-UV (A = 0.1-100 /Ltm) spectra of bright Seyfert 
galaxies. In reality, as we have shown in this paper, the 
decreasing trend of the MIR/bolometric luminosity ratio 
with increasing Lboi largely arises owing to the decreas- 
ing torus angular size, f^Torus, whereas the iDisk/iboi 
fraction remains nearly constant. 

7.3. Comparison with the UV-X-ray luminosity relation 

A number of studies have found that the ratio of 
near-UV (~ 25001) to soft X-ray (~ 2 keV) luminosi- 
ties in type 1 AGNs increases with luminosity (e.g., 
Vignali et ail [ 2003: St rateva et all 120051: iSteffen et alJ 
20061: lYoung et aJJl2010t see, h owever, I Yuan et al.lll998t 
Gaskell et al.ll20o4 iTang et al.ll2007t) . This suggests that 
the disk/corona luminosity ratio increases with luminos- 
ity, in apparent contradiction to our finding that this 
ratio is approximately constant over about two decades 
in luminosity. 

Part of the explanation may be that, although the 
near-UV flux might be a good proxy of the bolo- 
metric luminosity of the accretion disk in luminous 
quasars containing very ma ssive black holes, it might 
be a poor indicator (se e also iVasudevan fc Fabianl [20071 : 
iVasudevan et al.l I2009T ) in lower luminosity AGNs with 
less massive black holes (such as Seyfert galaxies), be- 
cause the maximum of their accretion disk emission is 
expected to be located in th e extreme-UV rather than 
in the near- UV ba nd (e.g., iShakura k, Sunvaevi 119731 : 
iHubenv et al.l |200TT ) . Mid-infrared observations, such as 
used in the present study, make it possible to disclose 
the true bolometric luminosity of the accretion disk by 
measuring the luminosity of the obscuring torus, which 
serves as a calorimeter of the power radiated by the cen- 
tral engine. 
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Furthermore, there is probably no discrepancy at all, 
because our study probes relatively low-luminosity AGNs 
compared to the quasars used in the L(2keV)/£(250(M) 
studies. Indeed, assuming again a T = 1.7 power-law 
spectrum with an exponential cutoff at E{ = 200 keV, the 
10 42 -10 445 erg s _1 luminosity range effectively probed 
by INTEGRAL in the 17-60 keV band corresponds to a 
range of 2 keV-monochromatic luminosities of 9 x 10 23 
to 3 x 1 26 erg s Hz" 1 . As can be seen, e.g., from 
Fig. 4 in iSteffen et all (|2006lK such AGNs are located 
in the low-luminosity part of the L(2500A)-L(2 keV) 
diagram, where the data are consistent with L(2keV) 
being proportional to L(2500A), whereas the evidence 
for a decreasing luminosity trend of i(2kcV)/i(2500A) 
comes from more luminous AGNs with L(2 keV) > a few 
xlO 26 erg s _1 Hz -1 . This possible change or gradual 
evolu tion of the trend is further suggested by our finding 
(see £16.51) that for typical quasars producing the CXB, 
with L(2 keV) ~ 2 x 10 26 erg s" 1 Hz" 1 , the L D i s k/i C orona 
ratio is about twice that for the INTEGRAL sample of 
(lower luminosity) AGNs. 

We conclude that although the reported behavior of 
the L(2keV)/L(2500A) ratio may be indicative of the 
iDisk/icorona ratio decreasing with luminosity in the 
most powerful quasars, this trend might be weak or ab- 
sent in less luminous AGNs. Clearly, further investiga- 
tion of this issue is required. 

7.4. Implications for SMBH radiative feedback 

The results of the present work have implications 
for the role of AGN feedback in the co-evolution of 
SMBHs and galaxies. The observed correlations between 
the masses of SMBHs and p arameters of their hos t el- 
liptical galaxies/bulges (e.g., iTremaine et al.l I2002D are 
possibly caused by coupling of energy released by the 
accreting SMBH to gas inside and around its host 
galaxy. This energy can be supplied in eit her mechan- 
ical ( e.g.. iChurazov et al.l [20 02; King 2003) or radiative 
(e.g., ICiotti fc Ostrikerll2001f ) form. The lat t er po ssibil- 
ity was considered in detail bv lSazonov et al.1 j2004) , who 
showed that radiation from typical quasars producing the 
bulk of the CXB is characterized by a Compton temper- 
atur<£3 of fcT c 2 keV. From th e INTEGRAL/ Spitzer 
study we now find, using eq. [1] in lSazonov et~atl ((2004) 
and taking into account the scatter in the LDisk/^Corona 
ratio, that for local AGNs, kT c vari es between » 2 and 
6 keV (note that the quasar SED in lSazonov et al.ir2004l 
extends to MeV energies, while here we do not take into 
account any radiation emitted above 300 keV). 

Thus, radiation from a SMBH can photoionize and 
Compton heat ambient interstellar gas above the virial 
temperatures of even giant elliptical galaxies. There- 
fore, AGN radiative heating can indeed play an impor- 
tant role in the co-evolution of gala xies and SMBHs, 
as has been suggested previously (|Ciotti fc Os trikcr 



2001 



iSazonov et al.ll2005l;lProga et al. 2008; N ovak et all 



20111 ). As discussed below, some studies suggest that 
the relative luminosity of the corona becomes small for 
SMBHs accreting near the critical Eddington rate, which 

10 The Compton temperature is the temperature of a gas inter- 
acting with a radition field at which there is no net energy exchange 
by Compton scattering between photons and electrons. 



would cause the Compton temperature to be relatively 
low for such actively growing black holes. Therefore, 
there is a need for a detailed study of the dependence 
of T c on the Eddington ratio. 

7.5. Dependence of the corona-disk relation on the 
Eddington ratio 

iVasudevan fc Fabian! (|2007l [20091 ). using UV and X- 
ray flux measurements of AGNs and estimates of their 
black hole masses, found that the disk/corona luminosity 
ratio is significantly larger for SMBHs accreting close to 
the Eddington limit, i.e., having Lboi/^Edd 1> relative 
to objects with iboi/^Edd "C 1. However, this result 
is associated with significant uncertainty (as emphasized 
by the authors), because the estimation of the intrinsic 
(unabsorbed) luminosity of an accretion disk is difficult 
even using UV data. 

It might be possible to reconcile the near constancy of 
the iDisk/icorona ratio found in our wo rk (a similar con- 
clusio n also follows from the study by IVasudeva n et al.l 
120101 who analyzed IRAS infrared data for a Swift/BAT 
sample of AGNs) with the s trong dependence of this ra - 
tio on Lhni/^F,dri inferred by rVasudevan fc Fabianl (|2007l ) 
and IVasudevan &: Fabianl (j2009) if we take into account 
the fact that most of the INTEGRAL ( and Swift) AGNs 
have relatively low Eddington ratios ()Horunzhev et aLl 
2012). On the other hand, a few of our objects, in partic- 
ular the narrow-line Seyfert galaxies, have iboi/^Edd ^ 1 
but nevertheless occupy the same locus on the £corona~ 
LDisk diagram as the other objects (Fig. [8]). This seems 
to argue against a strong dependence of iDisk/^Corona on 
the Eddington ratio. However, there is also a possiblity 
that the effective solid angle of the obscuring torus ex- 
plicitly depends on (decreases with) the Lboi/^Edd ratio, 
which would affect our estimates of the accretion disk 
luminosity based on the infrared luminosity of the torus. 

In future work, an infrared-X-ray cross-correlation 
study based on the INTEGRAL sample of local (rel- 
atively low iboi/^Edd ratio) Seyfert galaxies comple- 
mented by a representative sample of high Lboi/^Edd 
ratio quasars could help to clarify how the disk/corona 
luminosity ratio depends on the Eddington ratio. 
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APPENDIX 

A. IRS SPECTROSCOPY VS. IRAC PHOTOMETRY 

To verify our conclusions about host galaxy contami- 
nation of IRS spectra, we can compare 8 fim photometric 
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Figure 11. Starburst/AGN 15 fim flux ratio inferred from IRS 
spectra vs. extended/compact 8 fim flux ratio determined by IRAC 
photometry in 12" and 2.4" apertures (with estimated starlight 
contribution subtracted from the extended flux). Filled circles de- 
note objects from the clean sample (with AGN dominated spectra). 
Non-detections (< 2a) in either of the two quantities are shown as 
upper limits of 20% (shifted a little apart from each other in the 
lower left corner of the figure for better visibility). Empty squares 
denote the 6 objects with star burst dominated IRS spectra, for 
which the starburst/AGN 15 (im flux ratio is of the order of 1:1 
or larger. The dashed line indicates the ~ 3 : 1 extended/compact 
flux ratio expected for a sum of the adopted starburst template 
HBrand l ct al. 200(|) and a typical AGN spectrum in the case of an 
extended star formation region. 

measurements in small (2.4") and large (12") IRAC aper- 
tures. The result should depend on the extent of a star- 
burst. If a star formation region is more compact than 
the small (2.4") photometric aperture, there should be no 
significant associated extended 8 /im flux. In the oppo- 
site case of a starburst extending over > 10", i.e., beyond 
both the large IRAC aperture and the IRS extraction 
aperture, each 10% starburst contribution to the AGN 
spectrum at 15 /*m should add ~ 30% to the ratio of 8 fim 
fluxes in the 12" and 2.4" apertures, as follows from the 
comparison of our adopted starburst spectral template 
(jBrandl et al.ll2006l ) with a typical AGN spectrum (tak- 
ing into account that there are strong 7.7 /im and 8.6 /im 
PAH features falling into the 8 /im IRAC filter) . In addi- 
tion, there might be a non-negligible contribution of host 
galaxy starlight to the IRAC extended flux at 8 fim. In 
fact, extended stellar emission is clearly seen by IRAC 
in most of our sources, including "pure" AGNs, at wave- 
lengths A < 5 /im (see Fig.[T]). The starlight contribution 
to the extended flux at 8 /im can be roughly estimated 
as the extended 3.6 /im flux multiplied by (3.6/8) 2 (i.e., 
assuming a Rayleigh- Jeans spectrum). 

Figure [TT] demonstrates that almost all of the 
sources fulfil the above expectation, namely the ex- 
tended/compact 8 /im-flux ratio is less than ~ 3 times 
the starburst/AGN spectral flux ratio at 15 /im. The po- 
sitions of the sources in the IRS-IRAC diagram probably 
reflect different angular sizes of their starbursts: the fur- 
ther an object is located below the 3:1 limiting line the 



more compact is its star formation region. Two sources, 
1H 1934-063, and NGC 3081, lie somewhat above the 
3:1 boundary in the IRS-IRAC diagram. This appar- 
ent discrepancy probably arises because infrared spectral 
shapes of both starbursts and AGNs should in fact differ 
from one object to another, hence there may be a sig- 
nificant scatter around the ~ 3 : 1 relation between the 
starburst fractional contributions at 8 and 15 /im. We 
conclude that there is good overall agreement between 
signatures of host galaxy contamination found by IRS 
spectroscopy and by IRAC photometry for our objects. 

B. INVERSE CORRELATIONS 

The derived £15 (£hx) relation, eq. ((3]), makes it 
possible to predict the HX luminosity for a given MIR 
luminosity. For some applications, one may also be in- 
terested in knowing the distribution of £hx for a given 
Li5fj,m- We have therefore also computed an inverse lin- 
ear regression, i.e., log-Lux as a function of log £15 
for our clean sample excluding NGC 4395: 



i^HX.43 — (0.81±0.13)£ 



0.98±0.08 
15 ^m, 43 ■ 



(Bl) 



The rms scatter of £hx around the mean trend is 
0.39 dex. As can be seen in Fig. [2] this relation is differ- 
ent from the dependence of £15 u m o n £hx- 

The dependence given by eq. (|B1[) is expected to be 
affected by the Malmquist bias. Indeed, our AGN sam- 
ple is hard X-ray selected. Therefore, if for a given MIR 
luminosity £15 ^m, there is a range of possible HX lumi- 
nosities, Lhx, the INTEGRAL survey would find more 
objects toward the higher boundary of this range than 
toward its lower boundary because the more luminous 
objects can be detected from larger distances and hence 
from a larger volume of the Universe. For example, if for 
a fixed £15 Mm , £hx varies from object to object from L\ 
to £2 = 4£i, then the survey will find 4 3 / 2 = 8 times 
as many £2 sources as £1 sour ces. The intrins i c cor - 
relation can be found following Vikh lini n et al.l (j2009f ) 
by shifting the observed one (eq. |B1| ) by A log £jjx = 
-3/2 x In 10 x a 2 = -0.53, where a = 0.39 is the mea- 
sured rms scatter in log£nx- This results in 



£hx,43 (corrected) = (0.24 ± 0.04)£ ( j' 5 ! 



0.98±0.08 /nnl 

M m,43 ' KP 1 ) 

We can also compute the inverse relation between £Disk 
and £corona for our clean sample excluding NGC 4395: 



£corona,44 = (0.67 ± 0.06)££f™ 5 . 



(B3) 



It is also affected by the Malmquist bias. The intrin- 
sic correl atio n can be found by shifting the observed 
one (eq. |B3j ) by Alog£corona = — 3/2 x In 10 x cr 2 = 
—0.33 (where a = 0.31 is the measured rms scatter in 

fog £corona): 



^Corona 



44 (corrected) 



(0.31 ± Q.03)^ 8 ™ 5 . (B4) 



Although the inverse relations, eqs. (|B2[) and (|B4I) . 
are formally correct, they should be used with caution, 
because the implemented corrections for the Malmquist 
bias are comparable to or larger than the intrinsic scatter 
in the correlations (see Figs. [5] and [5]). The inverse rela- 
tions could be obtained more reliably using a MIR-flux 
selected sample of AGNs. 
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Figure 12. Illustration of the effect of source variability on the 
Lcorona - i'Disk correlation. Shown are the long-term light curves 
of NGC 4151 and Cen A in the 3-20 keV (RXTE/PCA, filled 
circles) and 17-60 keV (INTEGRAL /IBIS, open squares) energy 
bands. The sampling time is 1 and 3 days for RXTE and IN- 
TEGRAL, respectively. Also indicated are the "average" INTE- 
GRAL fluxes, /jjxi an d the periods used for their estimation (thick 
crosses), as well as the dates of the /i5 Mm measurements by Spitzer 
(arrows). Both /hx and /is Mm have been used in the MIR-HX 
cross-correlation analysis. Since the MIR data represent the ac- 
cretion disk emission averaged over the preceding period of a few 
years or longer, which does not coincide with the period over which 
the average HX flux was determined, it is clear that there should 
be significant associated scatter in the i/Corona^Disk correlation. 



c. 



EFFECT OF X-RAY VARIABILITY ON THE DERIVED 
CORRELATIONS 



The coronal emission detected by INTEGRAL is ex- 
pected to be substantially variable on time scales much 
shorter than the characteristic time scale of MIR varia- 
tions, and this should affect the observed £corona _ -E'Disk 
correlation. Indeed, Seyfert galaxies are known to be 
strongly variable in X-rays on time scales as short as 
minutes and more so on tim es scales of months and years 
(see, e.g,. IUttlev et al.l l2002). To roughly estimate the ef- 
fect of this variability on the ic rona _ ^Disk correlation, 
we have constructed long-term (1996-2009) light curves 
in the 3-20 keV and 17-60 keV energy bands of several 
bright, frequently observed AGNs from our sample using 
RXTE/PCA and INTEGRAL/IBIS data, respectively. 

Figure [T^] shows the light curves of NGC 4151 and 
Cen A, the two brightest AGNs in our sample. The 
(longer and better sampled) RXTE light curves demon- 
strate that the 3-20 keV fluxes of both sources varied 
by an order of magnitude over a time span of a decade. 
The same variability behavior is confirmed by the IN- 
TEGRAL data taken since 2002. The X-ray and hard 
X-ray variations are quantified in Figs.[T3"land ll41 which 
show the distribution of individual flux measurements, 
binned in log of flux. The observed flux distributions for 
NGC 4151 and Cen A can be well described as log-normal 
with a = 0.22-0.24 and 0.17-0.19 dex, respectively. This 



Figure 13. Binned distribution around the log-mean value of 
individual flux measurements for NGC 4151 by RXTE/PCA in 
the 3-20 keV energy band (upper panel) and by INTEGRAL /IBIS 
in the 17—60 keV band (lower panel). The solid line shows the log- 
normal distribution with a standard deviation equal to the actually 
measured rms scatter, 0.22 dex in the 3-20 keV band and 0.24 in 
the 17-60 keV band. 



Cen A 



10 - 




17 — 60 keV Flux, mCrab 



Figure 14. As Fig. 1131 but for Cen A. The rms scatter around 
the mean log F is 0.17 and 0.19 dex in the 3-20 keV and 17-60 keV 
energy bands, respectively. 

result is consistent with a number of previous studies that 
have demonstrated that X-ray flux variations in AGNs, 
similarly to X-ray b inaries, can be described by a l og- 
normal distribution (|Gaskellll200l lUttlev et al.H2005l ). 

The majority of AGNs in our sample were observed by 
INTEGRAL only occasionally, typically once a year for 
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a duration of a few days, hence the HX fluxes measured 
during 2002-2006, which are used in our cross-correlation 
analysis, in fact represent random snapshots of sources 
rather than their long-term averaged fluxes. Therefore, 
if NGC 4151 and Cen A are typical of the whole sam- 
ple, X-ray variability should induce a scatter of ~ 0.2- 
0.25 dex around the mean icorona _ iDisk trend. To fur- 
ther illustrate this point, HX fluxes averaged over a more 
recent period of INTEGRAL observations, from 2006- 
2009, give a factor of ~ 2 smaller Lcorona for NGC 4151 
and a factor of ~ 1.5 larger Lcorona for Cen A. 
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